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ABSTRACT 


This thesis describes the development of computer algonthms for experimental 
measurement control and subsequent signal processing of transient signatures to syn- 
thesize scattering impulse responses of scale model targets. The theories behind tran- 
sient scattering are considered in order to construct the algorithms. Up-to-date 
hardware and software technology are selectively 1mplemented to optimize the resultant 
signal-to-noise ratio. The detailed explanation of the hardware and software operation 
is provided. A noise model for the system is also discussed. Measurement and pro- 
gramming validations are described, where comparisons are made with numerical com- 
putations for transient scattering by selected canonical targets. Results are documented 


for futher research in the transient scattering problem. 
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I. INTRODUCTION 


A. BACKGROUND 

In circuit theory, it is well known that the characteristics of a linear time-invariant 
system can be obtained by exciting the circuit with an impulse function. The theory of 
transient scattering measurements 1s analogous to this method, and it also provides a 
unique parametric description of the scatterer via the complex pole locations in the 
Laplace domain transfer function. There are two approaches to obtaining the impulse 
response. One utilizes frequency domain methods and the other employs time domain 
techniques. The former has to sweep frequencies over the entire band but the latter uses 
direct pulse excitation. As a result, the latter can save time and expense. Therefore, 
direct transient scattering measurements are a viable alternative to the more conven- 
tional CW method for obtaining impulse responses of targets. 

In 1965 Kennaugh and Moffatt [Ref. 1] investigated the properties of the impulse 
response waveforms backscattered from metallic targets. The approximations for tran- 
sient scattering and impulse responses they made formed the basis for “ramp-response” 
imaging. 

In 1967 DeLorenzo [Ref. 2] presented a system that could be used to obtain an ap- 
proximate measurement of the impulse response of scattering bodies. His major devel- 
Opment was an indoor laboratory scattering range which facilitated the time domain 
analysis. He showed the agreement between his calculations and the measurement of the 
impulse response of a rectangular plate. 

Also in 1967, Murchison and Falk presented a hybrid system for measurement of 
impulse responses [Ref. 3]. This measurement system was separated into two frequency 
ranges. For frequencies below 4 Ghz, they used an impulse generator while for the fre- 
quencies above 4 Ghz, they used a swept frequency CW source. 

In 1971 Baum [Ref. 4] demonstrated theoretically that characteristic resonances can 
be obtained for scattering objects. His technique was called the Singularity Expansion 
Method (SEM). It was postulated that the resonances could be obtained from the 
transient scattering measurement for a complex object. 

In 1973 Tesche [Ref. 5] applied the SEM to to numerically evaluate complex reso- 


nances of simple metallic scattering objects. 


In 1975 Blaricum and Mittra [Ref. 6] proposed a technique for extracting the poles 
and residues of a system directly from its transient response. Their technique provided 
a number of numerical advantages. 

In 1980 the original Transient Electromagnetic Scattering Laboratory (TESL) be- 
came operational at the Naval Postgraduate School (NPS). Hammond [Ref. 7] devel- 
oped this laboratory for inverse scattering research based upon synthesized ramp 
responses. This laboratory utilzed an outside ground-plane range. This range was 
limited to symmetric scattering objects which could be “mirror imaged”. 

In 1983 the free-field TESL began operation at NPS [Ref. 8]. This laboratory is lo- 
cated in Spanagel 535. The TESL facilitated research efforts in high quality transient 
scattering measurements [Ref. 9] used for radar target identification and characterization 
using complex pole natural resonances [Ref. 10]. 

Prior to the current undertaking, McDaniel [Ref. 11] showed good agreement be- 
tween the computation and the measurement of the impulse responses for a sphere and 
a thin wire. These measurements and signal processing were performed using the 
Tektronix hardware set up, as originally employed in the 1980 TESL. In 1987 a rev- 
olutionary new Hewlett Packard digital processing oscilloscope became available. In 
addition, computational hardware had evolved dramatically since the early 1980's. This 
thesis describes the implementation of this new technology to enhance the performance 


of transient scattering and signal processing at the NPS TESL. 


B. OVERVIEW OF THIS EFFORT 

The main objective of this thesis was to incorporate new measurement and compu- 
tational hardware and to design and implement new software to automate high-quality 
transient scattering measurements. In this thesis, much work has been done in the de- 
velopment of the transient scattering work station which has the ability to measure and 
synthesize the scattering impulse response of scale model targets. The major hardware 
includes a microcomputer, a digital programable oscilloscope (DPO) and a shielded 
anechoic chamber. Automatic transient scattering measurements have been 1mple- 
mented by utilizing an [EEE-bus controller. 

With the new HP54120T DPO, the impulse responses of scale models are synthe- 
sized at a high signal-to-noise ratio (SNR). The acquisition and signal processing algo- 
rithms run much faster on the IBM PC-AT vis-a-vis the Tektronix interpreted BASIC 
computer previously used. The programs now use QUICK BASIC and FORTRAN 


compilers. Validations of software were made by way of comparing experimental 


measurements with computations for simple canonical targets. Various scale models 
were measured and documented in detail. Furthermore, a noise calculation was provided 
for further research. 

This thesis is broken into six chapters. Chapter II describes system facilities and 
hardware. This includes laboratory description, anechoic chamber, HP54120T digital 
processing oscilloscope and the GURU IIT hardware interface. 

Chapter III presents the theory of transient scattering measurements. This 1s com- 
posed of the free field range system representation, mathematical model and transient 
response solution. 

Chapter IV deals with acquisition algorithms and signal processing. This chapter 
covers operational description, GURU II software interface and the deconvolution al- 
gorithm. In addition, some selected outputs are illustrated. 

Chapter V emphasizes the experimental results. Two types of targets are presented. 
One is a canonical target and the other is a scale model target. The noise estimation of 
the system 1s also discussed in this chapter. 


Chapter VI summarizes all work and provides some future considerations. 


Il. SYSTEM FACILITIES AND HARDWARE 


A. LABORATORY DESCRIPTION 
Figure I on page 5 shows the general layout of the scattering range. It 1s composed 

of 

e Digital Programmable Oscilloscope (DPO) 

e Step pulse generator 

e Low loss cables 

e Broad band solid state power amplifier 

e Broad band transmitting and receiving antennas 

e Anechoic chamber 


¢ Personal computer (PC). 


Note that the DPO includes the sampling and pulse generator unit. The DPO 1s 
used as both the master pulse source and as the receiver. The antenna connection cables 
were obtained from the DPO accessory kit for use in low loss network measurements. 
They are flexible for hard-to-access areas. The characteristics are shown in Table 21 and 
Table 22 on page 78 in Appendix C. 

Transmitting and receiving antennas are of the double-ridged horn type. Their fre- 
quency range is 1-12.4 GHz, with less than I1dB of ripple in the boresight gain over this 
multi-octive band. The Avantek solid state power amplifier has an advertised 3dB 
bandwidth of 2 to 6 GHz but, by our observations, has a usable bandwidth of | to 7 
GHz. A I5V supply is needed for the amplifier. 

The anechoic chamber serves two purposes. First, it is used to isolate the target 
signal from the electromagnetic noise and interference outside of the chamber. Secondly, 
the absorbing material on the walls, floor and ceiling of the chamber are designed to 
provide minimal reflections, thus approximating unbounded space. The computer, 
which serves both acquisition and signal processing rolls, is an IBM-AT with 640 K 
RAM, a 20 megabyte hard disk, and a GPIB-PC2A interface card. The GPIB-PC2A 
was added to the PC in order to have the ability to talk or listen to the HP54120T DPO. 
The PC was set to run at 9 MHz. It has two floppy disk drives, and two output ports. 
One 1s a parallel port for the printer and the other 1s a serial port for the plotter. The 
color graphic board is NEC GBI. The displav is NEC EGA color monitor. 
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General layout of the scattering range 


Figure 1. 


B. ANECHOIC CHAMBER 

The anechoic chamber is located in room 535 of Spanagel Hall. This chamber in- 
corporates a metallically shielded anechoic enclosure to minimize outside EM interfer- 
ence. Chamber dimensions are 20 ft. long with a 10 ft. width and 10 ft. height. Because 
targets are suspended within the chamber, there is no requirement that they have mirror 
symmetry, as is the case for ground plane based ranges. In the measurements that were 
taken for this thesis, the target 1s mounted on a thin styrofoam column at a distance of 
8.5 feet from the transmitting and receiving antennas. Both receiving and transmitting 
antennas are placed close together to closely approximate monostatic measurements. 
The chamber interior 1s covered by a special carbon impregnated foam absorbing mate- 
rial. The material is cut into longitudinal wedges along the side-walls, floor and ceiling 
to act as conduits of EM energy towards the back wall. The back wall is covered by 18 
inch long pyramids to act as slow absorbing transitions for incident fields. Further in- 
formation about this free-field TESL can be found from Refs. 9, 11 and 12 


C. HP54120T DIGITAL PROCESSING OSCILLOSCOPE 

The Hewlett Packard 54120T is a fully programmable 20 GHz, four channel, 
digitizing oscilloscope with a nine inch color display. It is capable of automated meas- 
urements, digital storage, and TDR measurements. This DPO has only been in pro- 
duction since mid-1987 and represents the state-of-the-art in sampling technology. 
Previous to this thesis effort, a 6 year old 12 GHz DPO from Tektronix was used to 
perform transient scattering measurements in the TESL. 

The HP-DPO 1s comprised of two major components. These are the HP54120A 
Digitizing Oscilloscope Mainframe and the HP 54121A Four Channel Test Set. It also 
contains extensive self-tests to ensure proper functioning. Both HP 54120A and HP 
S4121A were calibrated together as a system called the HP 541201. The DPO can be 
remotely programmed by an HP-IB (IEEE-488) controller interface. The HP-IB com- 
plies with the IEEE 488.2 standard. The essential characteristics of the HP 54120T are 
summarized in Table 1 and Table 2 on page 7. 

It should be noted that the HP 54121A 4 Channel Test Set is sensitive to 
electrostatic discharge (ESD). Precautions should be taken to keep ESD from damaging 
the equipment. Details can be found in Ref. 13. 

Waveforms that are stored in the DPO memory can be sent to either an HP Thinkjet 
printer or an HP Plotter, both of which are accessed via the HP-IB bus. There are two 


ways to operate the instrument. One is manual operation via the front panel and the 


Table 1. HP54120T CHANNEL CHARACTERISTICS 


Bandwidth 20/12.4 GHz 
Number of channels 












Table 2. HP54120T TDR AND TIMEBASE CHARACTERISTICS 


less than 45 psec at 12.4 GHz bandwidth in average 
mode 


less than 10 psec + 0.1 % of reading 











other is to operate by sending HP-IB bus commands via the computer. A major effort 
in this thesis was developing programs to provide easy operator control and signal ac- 
quisition of the DPO using the computer. These algorithms are explained in detail in 
chapter IV. 

Ensemble waveform averaging has been implemented in acquiring scattering signa- 
tures. The averaging reduces noise by 1/,/N , where N is the number of averages. The 
minimum discernible signal level in the digital processing oscilloscope system is approx- 
imately 35uV. Typical noise amplitudes for the DPO sampling unit are shown in Table 
3 on page 8. 

Data display resolution is 1024 points horizontally x 256 points vertically (for dis- 


play times of at least 200 ps/div). The display can be selected in persistence or average 


Table 3. TYPICAL NOISE 












mode and the number of averages can be specified as a power of 2, up to 2048. The data 


Output transfer rate and the data record lengths are shown in the following table. 


Table 4. DATA OUTPUT TRANSFER RATE AND DATA RECORD LENGTHS 


10 psec/div < time/div < 20 psec/div 
20 psec div = time din = 20m sec ain 
50 psec/div < time/div < 200 psec/div 
200 paeerdiv = timed < 1 500m 






















For the time domain reflectometer (TDR) system, the risetime 1s an adjustable value 
based upon the timebase setting. The minimum is 10 ps or .08x time/div, whichever 1s 
greater. The maximum 1s 5 x time’div. Usually channel | 1s used for reflection meas- 
urements, however, 1n this thesis, it was used as the source generator. 

The digitizer uses a 12-bit successive approximation A/D converter. The full-scale 
range of the A’D is 640 mV. The least significant bit (LSB) of the A’D converter equals 
250 uV. This gives one part in 2560, or more than 11 bits of resolution. Averaging can 
stretch out the resolution to 32 uV. This improves the resolution to around 14 bits. 

The signal is sampled and digitized at a rate determined by the trigger rate, repe- 
tition rate, time base range, display mode, and number of channels turned on. Further 
information concerning the characteristics of the HP 54120T can be found from Ref. 
14. 


D. GURU II HARDWARE 
The Textronix supplhed General-purpose interface bus User’s Resource Utility 


(GURU II) 1s a combination hardware and software package that allows the use of any 


of the IBM personal computers (PC), and some PC compatibles, to control GPIB 
programable instruments. In this section, the focus is on the hardware. The software 
is explored in the next chapter. 

The hardware consists of the PC2A GPIB interface card that plugs into the PC main 
board, and the GPIB cable. The GPIB PC2A interface card is an 8-bit half-size board 
which performs as a talker, listener and controller. It also can be implemented with the 
full range of talker, listener, serial poll, service request, and remote programming func- 
tions for as many as 14 devices on the bus. It can perform as a complete controller. 
Using the PC2A, data transfer between the system memory and the GPIB are performed 
at rates of more than 300 kilobytes per second. The PC2A is compatible with all revision 
levels of the IEEE-488 standard including the HP-IB. 

The GPIB-PC2A interface card consists of these major sections: 
e address decoding 
e buffering and data routing 
¢ interrupt arbitration 
e direct memory access (DMA) arbitation. 
e configuration switches and jumpers 
e GPIB-adapter-TLC (talker/listener/controller) 
e time of day clock with battery back-up. 


The address decoding monitors the address lines to recognize when the GPIB-PC 
1'O address 1s present on the computer [/O channel and enables read and write access 
to the GPIB adapter register. The buffering and data routing handles data transfer be- 
tween the IBM PC-AT I/O channel and the GPIB adapter through a bidirectional 
internal data bus. 

Interrupt arbitation recognizes when interrupts have been enabled or disabled and 
passes or inhibits them accordingly. DMA arbitation recognizes when DMA operations 
are enabled or disabled, and when the last transfer has taken place. It also routes the 
DMA request and acknowledges signals to the selected DMA channel. The DMA 1s 
used to transfer data directly between a peripheral and memory without CPU. 

The IEEE-bus controller forms an interface between the computer and the 
IEEE-488 bus. Conceptually, the IEEE-488 bus behaves in a way very similar to the 
system bus. Hewlett-Packard originally devised it [Ref. 15] and intended it to link to- 
gether programmable instruments in a laboratory or industrial environment. By con- 


trolling test equipment and measuring devices from the IEEE-488 bus, implementing an 


automatic testing station 1s possible. A system under test is connected to the test 
equipment and measuring devices. The computer configures all equipment via the 
JEEE-488 bus and then reads the test results from the same bus. The GPIB has the 
following limitations: 

e Half of the devices attached to the GPIB must have power on. 

e The total length of the bus cable must not exceed 20 Meters. 


¢ No more than 15 devices may be attached to one bus. 


More detailed information can be found in Ref. 16. 
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I. THEORY OF TRANSIENT SCATTERING MEASUREMENTS 


A. FREE-FIELD RANGE SYSTEM REPRESENTATION 

The physical model of the system free field scattering range was first presented by 
Morgan in Ref. 12. The system block diagram is shown in Figure 2 on page 12. Im- 
portant components and parameters of the system are summarized in the frequency do- 


main representation in the following table. 


Table 5. FREQUENCY DOMAIN FUNCTIONS 


The multiple scattering between calibration sphere and 


absorber 














The transient pulse source 1s the step generator in the TDR system. This generator 
not onlv supplies the low-level transmitted pulse, before amplification, but also acts to 
coherently trigger the sampling circuits of the DPO receiver. Jt thus is akin to the 
“master oscillator” in a coherent pulsed CW radar system. The risetime of the step 
generator is approximately 25 psec. This step pulse is first amplified and filtered by a 
broad-band GaAs FET amplifier. 

The amplified pulse, which no longer resembles a step waveform, is then fed through 
the transmitting antenna which radiates into the anechoic chamber. The radiated field 
emanating from the antenna is approximately an angular weighted spherical wave. At 
the target position the incident field is a quasi-plane wave. Because of the close prox- 
imity of the transmitting and receiving antennas, some strong coupling of energy is in- 


duced in the receiving antenna. 
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This undesirable sigial can be seen clearly if the time delay setting on the DPO 1s ad- 
justed to begin many nanoseconds before the target response time-window. Therefore, 
the time setting delay must be checked when the target is present. | 

The initial calibration of the beginning of the target time-window was performed by 
using a rectangular copper plate, | ft wide and 1.5 ft high, as the test scatterer. The back 
scattering was observed with the copper plate turned broadside to the centroid of the 
horn positions and placed where the target would be mounted. The copper plate scat- 
tering signature was so large that it could be clearly discerned. The results of this 
measurement indicated a DPO ume delay setting of 38 nsec was needed. 

Multiple scattering between the target and the chamber absorber 1s represented by 
the two way interaction arrow connecting Hf) to H,(/). The total noise of the svstem 
N(f), is due to both thermal emission from the chamber and receiver noise in the sam- 
pling front end of the DPO. These sources are combined together and denoted as N,(/) 
for the n-th measurement that 1s made. To understand the theory of transient scattering 


measurements adequately, the mathematical model will be presented in the next section. 


B. MATHEMATICAL MODEL 

The five key parameters in Table 6 on page 14 play a major role in both the physical 
meaning and mathematical sense of the acquisition and signal processing algorithms. 
These five parameters can be related to the systein representation, described in the pre- 


vious chapter, by the following three equations: 


Vin = HA) HN AO LL + HN + Hs) + Hs] + M0 (3.1) 
YN = HN) He) X(N CL) + UN) + N40 (3.2) 
V3) = Hy) Ha X10) [Ha + He) + Hs) + Hs] + Ns. (3.3) 


These three equations represent the three respective measurements that are neces- 
sary for processing the scattering signature of a given target. As shown below, these 
measurements are for a calibration sphere, whose measured scattering can be compared 
to an accurate calculation; the “background” which denotes no target present, and the 
actual target scattering measurement. 

Equation 3.] states that the response to the amplified pulse generator, X(f), due to 
the calibration sphere and the chamber is equal to the product of the pulse excitation, 


transnutting antenna transfer function /7,(/), receiving antenna transfer function /7,(/) 


3 


Table 6. MATHEMATICAL NOTATION 


The subtracted target waveform 


times the sum of the antenna coupling H,{/) plus the calibration sphere H7S(f) and the 











multipath multiple scattering H¢(f). Added to this product is the thermal noise N,(/). 

Equation 3.2 represents the background measurement, with no target present. The 
svstem response is sinular to that previously considered, but with the absence of the 
transfer function for the calibration sphere. Equation 3.3 expresses the target measure- 
ment in the {frequency domain. The intepretation follows the same pattern as equation 
3.1 except the target transfer function, //{U/), appears instead of that for the calibration 
sphere. The multipath multiple scattering i equation 3.3 also changes and is denoted 
by Mgcf/). 

The three equations as stated above are three crucial measurements in the transient 
scattering problem. The backscatter measurement with the canonical (calibration) ob- 
ject },(Y/), measurements without target (background) Y,(/) and measurements with the 
desired target },(f) are obtained directly using the time domain acquisition program in 
Appendix B. This program was written in BASIC as part of this thesis effort. To speed 
up the execution of the program by a factor of about fifty times vis-a-vis direct inter- 
pretation, it has been compiled using the Microsoft Quick Basic Compiler. 

After acquiring accurate data, the next step is to process the data and extract the 


impulse response. This is the topic of the next section. 


C. TRANSIENT RESPONSE SOLUTION 

The three measurements that were just considered form the basis for estimating the 
“smoothed” impulse response of the target. The band-linuted nature of the measure- 
ments does not allow synthesis of the true infinite bandwidth impulse response which 1s 
due to an incident delta function plane wave. Instead, the “smoothed” term means that 
we can estimate the transient target response due to some predefined incident signal (a 
double-Gaussian shaped waveform is used here) whose significant spectral content does 


not exceed the bandwidth of the measurements. Aside from the currently used solid- 


state power amplifier, the bandpass range of the DPO-antenna system is about | GI1z 
to 12 GHz. The spectral bottleneck is due to the GaAS FET amplifier, whose passband 
is about 1 GlIz to 7 GHz. | 

In order to synthesize the smoothed transfer function of the target HZ(f), the first 
step is to subtract the “clutter” portion of the responses in equations 3.1 and 3.3 which 


is embodicd in equation 3.2. This subtraction results in equations 3.4 and 3.5 below: 


i= Y= Y2(0 (3.4a) 
Ya) = Hi) He) XI) [Hs + Us] + Nal) (3.46) 
Ni) = N,() — Ni) (3.4c) 
Y(N=%3N-¥2/ (3.5a) 
¥s() = Hy) Ua XD [sO + Msc] + Ns) (3.50) 
Ns(f) = Nx) — NV). (3.5c) 


Equations 3.4 and 3.5 are the very critical steps in obtaining the smoothed transient 
response of the target. The critical nature of the subtractions has to do with coherency 
in ume of the subtracted clutter. The calibration sphere measurement }/(/) had the 
background measurement }¥,(/) subtracted in order to eliminate the clutter. The clutter 
signal is that which 1s present aside from the target and is composed of two niain con- 
stituents. One is the direct coupling of the receiving and transnutting antennas, //,(/), 
and the other is direct scattering (as opposed to multipath) from the absorber in the 
chamber, //-(f). Since the measurements are made directly in the time domain, the sub- 
tractions are done in time domain. Ina simular fashion, equation 3.5 can be obtained. 

The results, },(2) and ¥;,(1), are then transformed into frequency domain by discrete 
fast Fouricr transforin( FFT ). This leads to the next step, which is to form an optimal 


deconvolution estimator for target response given by 


wise) 


Xo) He(f) = 
ie QP +c 


(3.6) 


where },(f) is the complex conjugate of Y,(/), Y;-(/) is the computed calibration sphere 


and C is the sinoothing paraineter. 


Ysclf) = Xol(f) Hs(N 3.7) 

X(t) — FFT > Xf) (3289 
2 i 

xo(t) = Ce Ot Cy eT A (3.9) 


where, for a specified o, and o, in this “Double-Gaussian” ( DG ) pulse, the C, and C, 


are chosen such that 


Xollo) = | (3.10) 
and 

+00 
[__,s(ndr = 0 (3.11) 


Each of the two Gaussian waveform in equation 3.9 will be specified here by a “10% 
pulse width”, defined by Az, such that 


CWA» 


aaa (3al2@) 
which gives 


4 9.21 
In(10) = — (3.125) 
Gay (A ty)" 


26s 








The estimator in equation 3.6 1s due to Riad [Ref. 17] and results from an optimal 
least-square fidelity criterion for deconvolved target response. The result of equation 
3.6, X,(f) Hi , represents the frequency domain scattering response of the measured tar- 
get to an incident plane wave whose Fourier transform is X,(f). This incident waveform 
is specified by numerically computing the scattering response of the calibration sphere 
to it. This numerically generated response of the calibration sphere is embodied in 
Y;-(f) , as shown in equation 3.7. It is obtained from the SPRSCT program which uses 
the magnitude and phase data from the MIE series program. Both of these programs 
were written by Morgan. 

Both subtraction and optimal deconvolution algorithms are implemented in the 
Deconvolution Program in Appendix A. The procedure to operate this program will be 
explained in chapter IV. 


IV. ACQUISITION AND SIGNAL PROCESSING ALGORITHMS 


A. GURU IIT INTERFACE 
1. GPIB Principal 
This section emphasizes the GURU II software interface. In order to under- 
stand the Acquisition program operation, the GPIB principle has to be understood. The 
General Purpose Interface Bus uses a twenty-four conductor cable. Sixteen of these are 
used to send the signal. Data its sent back and fourth between the listener and the talker 
in a bi-directional fashion. This means that both transmitted and received data travel 
on the same line. Eight lines are used for 8 bits of parallel data transfer. A byte 1s sent 
by serial mode but each bit 1s sent in parallel mode. Three lines act as the coodinator 
between the talker and the listener. These lines are assigned their duty as follows: 
e no data available line 
e data not ready line 


e¢ data not accepted line. 


The talker on the bus controls the data available on the line and waits for the 
listener's response on the “not ready for data” line or “data not accepted” line. The re- 
maining five lines are the attention line, the interface clear line, the remote enable line, 
the service request line and the end or identify line. When the attention line is active, 
the talker and the listener will be specified. The controller has full access to five of these 
lines and has the capability to directly address other devices to command them to talk 
or listen. The interface clear line is used to set the turn-on condition. When the remote 
enable line is asserted, all devices on the bus are under computer control and nothing 
external can interfere. The service request line is used by a device to indicate a require- 
ment for attention and to request an interruption of the current sequence of events. 
There are two methods for the controller to determine the device and the message. One 
is the serial poll while the other is the parallel poll. The polls obtain numerical codes 
which show a specific event. The serial poll is a mechanism that the controller uses to 
obtain status from individual devices. As each device is polled, it returns a status byte 
which controller decodes to determine the device requesting service. The parallel poll is 


a mechanism for accepting and decoding the general status of as many as eight devices 


Uy 


simultaneously. The end or identify line is used by the talker to indicate the end of the 
message. 
2. Basic Software Requirements 
There are 3 files which the GURU requires in order to run the application pro- 

gram. These files are 

¢ GPIB.COM 

© IBCONF.EXE 

e CBIB.OBJ. 


The primarv GPIB device drivers are read in via the CONFIG.SYS file at boot-up time. 


This file is in the root directory and must include the following DOS statement: 
¢ DEVICE=GPIB.COM. 


This causes the GPIB driver in GPIB.COM to be loaded when the system is booted. 

The IBCONF.EXE is the executable program to configure the GPIB.COM file 
containing the device drivers for each instrument. HPDEV is the name of the device 
driver for all Hewlet-Packard instruments, including the HP 54120T, whose decimal pri- 
mary address is 17. The IBCONF.EXE (Interface Bus Configure) program allows a 
match between the softwere and the hardware. The program prompts with menus of 
details about the hardware parameters associated with the interface card. This device 
for Hewlett-Packard is already configured in the system. When the IBCONF.EXE is 
run, the device drivers are mapped for each instrument in the display. By selecting the 
HPDEV, and editing it, the characteristics of the HPDEV show on the screen. This al- 
lows the operator to change the characteristics to be compatible with the hardware. 
After reconfiguration of the GPIB.COM file, the operator must return to the operating 
system and then re-boot the PC for the new parameters to be in effect. 

CBIB.OBJ is the object file which is linked with any compiled BASIC program 
which needs to interface with the GPIB.COM device drivers. In our case, the Acquisi- 
tion Program, written in GW-BASIC, 1s first compiled using the Microsoft QUICK 
BASIC compiler. This process creates an object (.OBJ) file. The CBIB.OBJ object file 
is then linked with the Acquisition object file. The final result is an executable (.EXE) 
file. 


B. ACQUISITION ALGORITHM 
The Acquisition program was written in Microsoft GW-BASIC. The hardware is 
interfaced with this program by way of the CBIB object file. This file comes with the 
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GURU II software. After compiling the source code, the object file of the Acquisition 
has to be linked to the CBIB object file in order to complete the creation of the Acqui- 
sition executable file. To run the program, the Brun3087 execution file is needed in the 
same DOS Directory. This file is supplied with Microsoft QUICK BASIC. The Ac- 
quisition program uses a friendly and interactive menu format. The program begins with 
a displaved welcome message, including the time and date. This screen also provides a 
brief procedure on how to measure the transient scattering response. The screen gives 
operator the choice of using the system default or the last response. The system default 
was designed for first time operators. The last response was designed for experienced 
Operators. A main menu 1s displayed on the screen after the operator decides whether 
to use the system default or the last response. 

The main menu contains many importance parameters. A first item is the number 
of data points in the acquired time series. This number has to be a power of 2, but not 
more than 1024. This limitation could be modified in the program by defining a larger 
array size up to 2048 points. The second item is the number of sub-averages. This pa- 
rameter determines the number of values to average for each time point when the 
HP54120T is in the average mode. This parameter has to be an integer from | to 2048, 
also in a power of 2. The third item is the number of data blocks. This determines the 
number of sets of data in one record. The fourth and the fifth items are a user supplier 
identification for the target and the date of measurement. The sixth item tells the oper- 
ator about the typical sampling time interval. The seventh parameter is the type of the 
waveform being acquired on the measurement being made. As discussed in the previous 
chapter, there are three scattering waveforms which need to be acquired for each target 
being considered: calibration (sphere), background and actual target. The eighth item 
tells the operator the time window in nanoseconds. The ninth, tenth and eleventh items 
are the maximum vertical setting, the data file name output and the auto time setting. 
At the bottom of the menu, the operator is asked if he wants to change the items. 

The second menu appears on the screen when the operator decides to change some 
item on the main menu. This “overlay” menu contains all changeable parameters and 
allows a change to these parameters interactively. When the changes are completed, the 
operator can go back to the main menu by typing the letter “R”. This brings the oper- 
ator to the main menu and displays the new item on the menu. When everything is 
ready, the operator selects “no change” to begin to automate the DPO. The program 


starts by initializing the DPO and setting all specified parameters as directed by the main 
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menu. The program stops when the DPO has finished the setting operation. The pro- 
gram then asks for the bandwidth of the measurement. There are two choices: 20 Ghz 
or 12.4 Ghz. To protect from errors in the initial setting, the program stops for the final 
check. As soon as the operator initiates the acqisition cycle, the beginning time shows 
on the screen. For a typical acquisition composed of 2048 averages and | data block, 
it takes about 7 minutes to finish the operation. In the meantime, the computer waits 
for the DPO to provide the raw data. When the DPO finishes this operation, the ending 
time is displayed and the PC opens a temporary hard-disk file for storing the raw data. 
Before the operation is ended, the header which is to be placed on the data file is dis- 
plaved on the screen. The operator can then observe any error in the data set, should 
there be any. 

The header 1s transferred as a separate file. The header provides important charac- 
teristics of the data, such as format type, number of points, number of averages, x in- 
crement, X origin, x reference, y increment, y origin, y reference and range setting. All 
of these parameters allow the program to recover meaningful data from the raw data. 
The remainder of the program constructs the standard data file structure and the hard 
copy output. 

At the end of the program, there is a pause with a query regarding the operator's 
intention to continue with an additional measurement. To complete the acquisition cy- 
cle for one target, the three measurements previously discussed have to be made in order. 
These three files will then be supplied to the deconvolution program in the next proc- 
essing step. For acquiring more than three measurements, the background must be 
measured for each target so it can be subtracted from the respective target measurement. 
The background measurement occurs either immediately before or after the target 
measurement. This reduces errors induced by long term nonstationary drifts in the sys- 
tem. The need for such temporal coherence of the target-background measurements has 
been verified experimentally during the course of this work. Typical outputs are as 
Shown in Figure 3 on page 21. 


C. DECONVOLUTION ALGORITHM 

The deconvolution program was written in Ryan-McFarland (RM) FORTRAN in 
order to improve the run time. The source code is standard ANSI 77 FORTRAN 50 it 
may be recompiled under virtually any PC FORTRAN compiler such as Microsoft or 
Lahey. This program was originally written in Tektronix BASIC by Morgan [Ref 11] 


and was restricted to use on a Tektronix 4052 microcomputer. 
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DECONVOLUTION INPUT STEPS 
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Figure 3. Deconvolution input waveform 


The use of a standard FORTRAN-77 compiler generalizes the program so it can be 
run on any IBM PC compatible. This algorithm follows the mathematic deconvolution 
as described in chapter III. The program has graphic capabilities which are supplied by 
GRAFMATIC software. The procedure to compile the program is described in the RM 
FORTRAN manual. This program must be linked to the GRAFMATIC library in or- 
der to run properly. 

The program begins by displaying the standard input data file. This data file format 
has been standardized for use with all of the software support systems. The header of 


each data file contains the important information such as identification, type of meas- 
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urement, date of measurenient, time of experiment, the maximum time scale, data block 
and number of points. Types of measurement are as explained in the previous section. 
The first data file needed for the deconvolution program is that for the computed seat- 
tering by the calibration sphere target. This computation data is obtained from the MIE 
series program and SPRCT program, written by Morgan, as described in the previous 
chapter. The input data required for these programms are summarized in Table 7 on page 
Zoe 

The sphere scattering computation 1s transformed into the frequency domain and 
displaved on the screen in order to observe the bandwidth. The bandwidth of the com- 
putation has to be approximately the same as the estimated bandwidth of the measure- 
ment. If the bandwidth is not appropriate, the operator should to terminate the 
program via a Ctl-C and run the MIE series followed by SPRCT program to provide 
appropriate data. 

In the next step, the program asks for the calibration and background file names. 
Both data are read in and plotted on the same axis in order to observe if time slufting 
has occurred, and ifso, to be able to estimate how much. The operator is then allowed 
to shift one waveform either Ieft or right and then re-observe the overlayed waveforms. 
Once a satisfactory alignment is observed, the calibration waveform is subtracted from 
the background waveform and the difference waveform 1s transformed into the frequency 
domain. 

The program then stops and asks for the target file name. It then asks if the oper- 
ator desires to use the same background or another, depending on which background 
was closer to the target completion acquiring time. In a similar fashion, the target data 
and the background data are plotted on the same axis. If there is an observed drift, the 
program allows the operator to shift the target waveform to the left or to the nght until 
both target and background are lined up. Then the subtraction can be performed. The 
difference waveform 1s then transformed into the frequency domain. Typical scattering 
data being processed by deconvolution program 1s shown in Figure 3 on page 21. The 
overlay and shift is shown in Figure 4 on page 23 and the subtracted waveform 1s shown 
in Figure 5 on page 24. The target for this example was a 10 centimeter long “thick” 
copper wire, and the calibration target was a standard 0.1025 meter radius aluminum 
sphere. After the subtraction, the difference calibration waveform and the difference 


target waveform were tranformed into the frequency domain. 


paps 


TARGET AND BACKGROUND MEASUREMENT COMPARISON 





aP 
O 
S ; 
© 
5 
oN TARGET 
fr eS everett BACKGROUND 
= oO 
a 
z 
& 
> 
=o 
Lo 
& 
E 
Ow 
Zz (e) 
w 
= 
{ 
- ‘ 
O 
~ 
G f 
0 § 10 15 20 


TIME IN NANO SECONDS 


Figure 4. Overlay 


Table 7. CALIBRATION SPHERE PARAMETERS 


Sphere radius in meters ORlO25 
Distance from sphere in meters 
Bistatic angle in degrees 


Scattering plane 
Time window in nanoseconds 
Nyquist frequency in Ghz 
No. frequency points 
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Figure 5. Subtraction waveform 


The three frequency domain waveforms were formed into Riad’s optimal deconvo- 
lution estimator given in equation 3.6. The smoothing constant “C” is selected by the 
operator. This constant is set to optimize the denominator threshold level in the fre- 
quency domain. If it is set too large, oversmoothing will result which reduces the spec- 
tral content in both the low and high frequency regimes. If the threshold is set to zero, 
a “naive” deconvolution results. This provides vastly amplified noise in the signal and 
typically looks like “garbage”. A typical value for Cis 0.1. Riad’s method yields an 
estimator for the scattered waveform of the target due to an incident plane wave signal 
having signal shape equal to that which illuminates the computed calibation sphere us- 


ing program SPRCT. This frequency domain estimator is transformed to the time do- 
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main and the result is displayed. At the end of the program, the operator may try 
another smoothing parameter. If the result satisfies the operator, it can be stored in the 
file whose file name was provided by the operator. The program also allows the proc- 
essing of another target directly without rereading the computed data and the calibration 


data. 


D. EXAMPLE OUTPUT 
The output of the deconvolution program 1s in both the time domain and frequency 
domains, but the following example shows only the time domain output. This output 


can be stored on a hard disk or a floppy disk. 
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Figure 6. Thick wire measurement 
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The impulse response shown in Figure 6 belonged to the signature of the thick copper 
wire having a 90 degree incident wave and a 90 degree reflected wave. The number of 
averages used in the measurement was 2048. The smoothing was C=0.1. There was 
no time drift observed. One data block was used. Other parameters can be found in the 


following table. 


Table 8. 90 DEGREE THICK WIRE PARAMETERS 


Wire length in meters 0.1 ae 
0.0026 
Incident angle in degrees 1/90 


DG pulse minimum 10% spread in nanoseconds 


DG pulse maximum 10% spread in nanoseconds 030 












0.28 
0.30 
[Time window in nanoseconds SSOSCSCSC~S~S~S~STOO 
[Scattering angle in degrees —~—~SCS~“~S*~“S*S*~*~*~S~*~dC 


Note that the minimum and the maximum 10% spreads determine the respective «, 
and a, in the DG waveform of equation 3.9. 

First, the frequency spectrum of the smoothed impulse response was calculated from 
the output and then transformed to the time domain by an FFT [Ref. 18]. The reason 
for calculating the frequency response first 1s that it 1s more convenient to work in the 
frequency domain when using Raid’s method to extract the smoothed impulse response 
using equation 3.6 

The dotted line in Figure 6 on page 25 is the actual measurement. Because the 
measurement was contaminated by noise, the result did not agree exactly. The contin- 
uous line represents the computed thick wire smoothed impulse response due to the 
Same excitation. This computation was performed directly in the time domain using an 
E-field integral equation. Measurements for other types of structures will be illustrated 


in the next chapter. 
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V. EXPERIMENTAL RESULTS 


A. CANONICAL MODEL SCATTERING 
In addition to the 10 cm long thick wire considered in Chapter 4, four additional 
canonical targets were Selected: 
¢ 10 centimeter thin wire (broadside and canted incidence angles) 
¢ §8.1 centimeter diameter sphere 
e 12.2 centimeter diameter sphere 


e 15.5 centimeter diameter sphere 


The Mie series program, MIE, was used to compute the frequency domain scatter- 
ing transfer function for the various metallic spheres. These transfer functions were then 
used as inputs to the SPRCT program to compute the transient scattering responses due 
to the specified Double-Gaussian (DG) pulse incident plane wave, as specified by 
equation 3.9. The thin wire computations were performed using the TDIG program. 
This program utilizes a numerically evolved time-domain electric field integral equation 
for the induced current distribution on the wire. The scattered field 1s then calculated 
by appropriate numerical space-time integrations using this induced current. This and 
the sphere scattering programs were written by Morgan. The measured and computed 
DG pulse scattering responses are overlayed in each plot to validate both the fidelity of 
the experimental procedure and the accuracy of the deconvolution algorithm. The 
characteristics of each canonical target scattering case are summarized in Tables 9-13 
on pages 28-33. A noteworthy item 1s that there is a slight bistatic angle (3 degrees) 
between incident and scattered aspects. This has been incorporated into the computa- 
tions. As can be readily observed, the agreement between measurements and computa- 
tions in Figures 7-11 on pages 28-32 is quite good in all cases. The additive “rippling” 
of the experimental results 1s due in part to the residual noise and clutter contributions. 
This latter pollutant includes uncancelled multipath interactions of the target and 
canonical sphere scatterers due to imperfect absorption by the anechoic chamber sur- 
faces. Additional differences in the comparisons are due to the approximations which 


are innate to the computations. 
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MEASUREMENT AND COMPUTATION COMPARISON 
12.2 CENTIMETER DIAMETER SPHERE 
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Figure 7. The 12.2 cm diameter sphere validation. 


Table 9. 12.2 CM SPHERE PARAMETERS 


Sphere radius in meters 0.061 
Distance from sphere in meters 
Bistatic angle in degrees + 














3 
Scattering plane 
Narrow 10% pulse width in nanoseconds 


Wide 10% pulse width in nanoseconds 
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MEASUREMENT AND COMPUTATION COMPARISON 
10 CENTIMETER THIN WIRE 


1.2 


0.8 
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AMPLITUDE IN VOLTS/METER 
-0.4 
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Figure 8. The 10 cm thin wire validation at 45 degrees 


Table 10. 45 DEGREE THIN WIRE PARAMETERS 
Wire length in meters 

Wire radius in meters 

Incident angle in degrees 

DG pulse minimum 10% spread in nanoseconds 

DG pulse maximum 10% spread in nanoseconds 


Scattering angle in degrees 


No. Wire segments 
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MEASUREMENT AND COMPUTATION COMPARISON 
10 CENTIMETER THIN WIRE 


0 0.5 1.0 
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Figure 9. The 10 cm thin wire validation at 90 degrees 


Table 11. 90 DEGREE THIN WIRE PARAMETERS 
Wire length in meters 

Wire radius in meters 

Incident angle in degrees 

DG pulse minimum 10% spread in nanoseconds 


DG pulse maximum 10% spread in nanoseconds’ - 


No. Wire segments 





Scattering angle in degrees 
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MEASUREMENT AND COMPUTATION COMPARISON 
8.1 CENTIMETER DIAMETER SPHERE 
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Figure 10. The 8.1 cm diameter sphere validation 


Table 12. 8.1 CM SPHERE PARAMETERS 


Sphere radius in meters 0.0405 
Distance from sphere in meters 
Bistatic angle in degrees 








200 
Nyquist frequency in Ghz 12.775 
Narrow 10% pulse width in nanoseconds 


Wide 10% pulse width in nanoseconds 
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MEASUREMENT AND COMPUTATION COMPARISON 
15.5 CENTIMETER DIAMETER SPHERE 
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Figure 11. The 15.5 cm diameter sphere validation 


B. SCALE MODEL SCATTERING 

The scattering responses of the wire and sphere targets due to the incident DG in- 
cident pulse can be numerically evaluated using straightforward algorithms. This is not 
the case for more complex scattering objects, such as a highly realistic model of an air- 
craft. Primitive computational algorithms do currently exist for such complex structures 
but the poor accuracy of their results would be insufficient to properly validate the 
measurements being made here. On the other hand, the very good comparisons shown 


in the previous section for both low-damped (wires) and high-damped (spheres) strongly 
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Table 13. 15.5 CM SPHERE PARAMETERS. 













support the conjecture that high quality measurements can be obtained for all metallic 
targets having comparable physical dimensions. 

After the validations were made using wires and spheres, extensive measurements 
were made of the transient backscattering from four different silver coated 1/72 scale 
plastic aircraft models. The full-size dimensions of each aircraft, whose identities will 
not be given here, are shown in Tables 14-17 on pages 34-37. 

Measurements were made for various incident aspects in the wing-fuselage plane on 
each model, using horizontal incident polarization: 0 degrees (nose-on); 30 degrees; 90 
degrees (broadside); and 180 degrees (tail-on). In addition, a broadside look-down as- 
pect was considered, with backscattering measured for incident linear polarizations both 
parallel to (“vertical wing”) and perpendicular to (“horizontal wing”) the fuselage. 
Samples of these numerous measurements are displayed in Figures 12-15 on pages 34-37. 

The complete library of waveforms is being used to support research in demon- 
Strating the accuracy and aspect independence of natural resonance radar target iden- 
tification. 

Four scale model aircraft were selected. They are as follows: 

e Scale model 1 
e Scale model 2 


e Scale model 3 


e Scale model 4. 


C. SYSTEM NOISE ESTIMATION 
Noise in the data acquired from the acquisition program was investigated in this 


effort in order to develop a mechanism to estimate the signal-to-noise ratio (SNR). The 
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Figure 12. Vertical and horizontal wing scale aircraft model | 


Table 14. DIMENSIONS OF SCALE AIRCRAFT MODEL 1 


Scale 
Length 36 feet 4 inches 


Wingspan 
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SCALE AIRCRAFT MODEL 1 


HEADING 000 DEGREE HEADING 030 DEGREE 


2 & 2 4 
TME IN NANO SECONDS TIME IN NANO SECONDS 
HEADING O80 DEGREE HEADING 180 DEGREE 


4 
TIME IN NANO SECONDS 


AMPLITUDE IN VOLTS/METER 
1.00 0 0.004 
AMPLITUDE IN VOLTS/METER 
—0.004 @ 0.004 


2 4 
TME IN NANO SECONTS 


Figure 13. Scale aircraft model | at 000, 030, 090, 180 degree aspect angle 


Table 15. DIMENSIONS OF SCALE AIRCRAFT MODEL 2 


7 
Overal nat 
Extended wingspan | 28 feet 7 inches 
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Figure 14. Scale aircraft model comparisons at 000 degree aspect angle 


Table 16. DIMENSIONS OF SCALE AIRCRAFT MODEL 3 
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HEADING O90 DEGREE 


SCALE AIRCRAFT MODEL 1 SCALE AIRCRAFT MODEL 2 
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Figure 15. Scale aircraft model comparisons 090 degree aspect angle 


Table 17. DIMENSIONS OF SCALE AIRCRAFT MODEL 4 


63.8 fee 


22.8 fet 
18.6 fee 
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SNR is a major determining factor for the ineasurement fidelity. Much of this SNR 
assessment involved the computing some statistics of the waveform. All of these statis- 
tical calculations were performed using GRAFSTAT with APL. 

In the acquisition step, the calibration, the background, the target and a second 
background were measured consecutively. The two backgrounds were measured for the 
purpose of noise estimation. In the following analysis, a 2048 ensemble-average mode 
was used to acquire the data. The acquisition of each waveform took about 9 minutes 
to complete. 

In a separate investigation, the stationarity of the system noise was investigated by 
acquiring seven consecutive background measurements. These backgrounds were then 
subtracted from each other to observe the noise occuring during the consecutive meas- 
urement. For example, the second background and third background were subtracted, 
as 1s shown in Figure 16 on page 39 along with various statistical parameters in Table 
18 on page 39 inentioned, the statistics of each subtracted waveform was calculated by 
the GRAIFSTAT software packgage. The variance of the noise can be calculated from 
the standard deviation and expected power can be approximated by this variance when 
the noise has a small mean. This approximation was reasonable because the squares of 
the mean for the actual measurement data are very small, ( see Figure 16 on page 39 ). 
The computed noise power is plotted versus time in Figure 17 on page 40. As can be 
seen, the measured noise power varied only slightly with time. A linear least-squares fit 


vields the equation 
Y = ~—93.042 + 0.035727% , (Say 


where Y is the expected noise power and X is the time in minutes. This model was 
based on only seven background measurements with only 2048 averages per measure- 
ment. For practical purposes the system noise can be assumed to be fixed in time, de- 
pending only upon the number of ensemble averages used to complete the waveform 
measurement. 

To calculate the SNR, the thick wire measurement was used. The signal plus noise 
waveform is displayed in Figure 18 on page 41 while the computed statistical parameters 
are given in Table 20 on page 41. 

This waveform was time windowed to eliminate the noise outside of the window 
before the statistical parameters were calculated. The signal plus noise power was cal- 


culated by way of the variance using the following equation. 
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EXAMPLE NOISE WAVEFORM 


SCATTER PLOT TABLE 


x 1 900K 
Y BGT 2R-BCT3 
2ND BACKGROUND MINUS SRO BACKGROUND SELECTION sALL 
ie x LABEL : TIME IN NANO SECONDS 
a Y LABEL : AMPLITUDE IN NANO SECONDS 
= NO. OF ELBENTS :812 
CORRELATION XY :70.34118 
1 RK CORRELATION 170.14703 Tec3.3568 
:40 
5 STD. DEVIATION :5.7904 
ERCENTILE  10.98230 
ERCENTILE : 4.9902 
z IAN 10 
of-- 75-PERCENTILE 118.01 
am ii GS-PERCENTILE :19.018 
= X MIN. :0 i 78 
& % MAX. :20 19.061 19.022 
. Y MZAN :~3.0211E71 
z STD. DEVIATION :0. 
6 B-PERCENTILE —1-0.000034289 
ILE :77.02E-6 
AN :-7.019E-8 
4 ILE :0.00002423 
x §S-PERCENTILE :0.000024231 
T Y MIN :“0.000088519 —0.000089519 —0.00008e519 
Y Wax :0.000086731 0.000055481 0.000085484 





Figure 16. Noise waveform statistics 


Table 18. STATISTICS OF NOISE WAVEFORM 


Mean -3.0211E-10 
0.000023345 


-7.019E-8 
Maximum 0.000066731 
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NOISE ANALYSIS 


SCATTER PLOT TABLE 





x > TIME 
Y : 10x10®POWERS 
SELECTION 2ALL 
X LABEL :>TIME IN MINUTE 
NOISE POWER OF SIX TRIAL no. oF Eten - 
° CORRELATION XY :0.60165 
RK CORRELATION :0.65775 T=#1.656 
X MEAN 731.5 
¢ STD. DEVIATION :16.857 
S-PERCENTILE 79 
25-PERCENTILE :18 
R : Y = -93.042 + 0.035727 x X re Beenie arene 
- ‘ Q@S-PERCENTILE 3:54 
X MIN. 79 18 27 
X MAX. 704 45 36 
Y MEAN 281.917 
STD. DEVIATION :0.999B6 
S-PERCENTILE 282.76 
; = _ = 25-PERCENTILE  :7-92.716 
TIME IN MINUTE MEDIAN 7792.3543 
- 7S-PERCENTILE :780.669 
9S-PERCENTILE :—90.669 
Y MIN 2782.76 -92.716 —92.636 
Y MAX :~90.669 —90.669 —92.049 


Figure 17. Noise power 


Table 19. STATISTICS OF NOISE MODEL 
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THICK WIRE SIGNAL PLUS NOISE 


SCATTER PLOT TASLE 
XXX 


x 
Y : TOTAL» ( TTHICK=TTHIBGT) 
SELECTION tALL 

X LABEL TIME IN NANO SECONDS 

Y LABEL sAMPLITUDE IN NANO SECONDS 
NO. OF ELEMENTS :512 

CORRELATION XY :0.0089981 

RK CORRELATION :0.0081907 T=0.13981 

X MEAN 210 

STD. DEVIATION :5.7904 

S-PERCENTILE  :0.98230 

25-PERCENTILE  :4.9902 

MEDIAN £10 

7S-PERCENTILE 315.01 

S5-PERCENTILE :19.018 

X MIN. :0 0.030139 D.078278 

X MAX. :20 18.961 19.822 

Y MEAN :76.2891E-8 

STD. DEVIATION :0.0003382 

B-PERCENTILE :-0. 00088722 

25-PERCENTILE  :0 

MEDIAN :0 

75-PERCENTILE  :0 

GS-PERCENTILE  :0.00059718 

Y MIN :-D.0014822 0.001308 -0.0012747 
Y MAX :0,0014128 0.001319 0.0012253 


— CC 


Figure 18. Signal plus noise waveform statistics 


Table 20. STATISTICS OF THICK WIRE SIGNAL PLUS NOISE 
WAVEFORM 


[Mean OOSOC—SCSCSCC“‘(s*‘O. 











4] 


E(X’) = Variance + E’(X) (5.2) 


where E(X) is mean of the noise and E(X?) is the mean square value of the noise (noise 
power). Note that the second term on the right side of equation 5.2 was close to Zero, 
sO it can be neglected. As a result, the expected power is equal to the variance. The 
signal plus noise power was equal to 0.0003382? or -69.41 dB. The noise expected power 
_ was calculated by equation 5.] and is given by -92.64 dB for the average time of 10 
minutes. Since the noise power was so small compared to the signal plus noise power, 
the subtraction between the signal plus noise and the noise was not necessary. The 
noise power when time windowed was reduced by a factor of 200/512, or 4.08 dB. The 
result for the SNR is 27.31 dB in the case of the thick wire measurement. Note that this 
SNR is about 10 dB better than that provided for the same target using the Tektronix 
equipment, as carried out by McDaniel [Ref. 11]. 
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VI. CONCLUSIONS 


A. SUMMARY 

This thesis describes the theory and implementation of a newly updated transient 
scattering measurement facility. An overview of the historical precedents and moti- 
vations for this development was given in Chapter J. This was followed in Chapter J] 
by a detailed description of the new equipment being used as well as the physical struc- 
ture of the anechoic chamber. 

Building upon this overview of the hardware aspects, the conceptual and math- 
ematical modeling of the scattering range 1s undertaken in Chapter IIJ. Using the 
physical model of the radiating and receiving antennas in the anechoic chamber, the 
various electromagnetic interactions between the scattering target and the chamber are 
represented by a comprehensive linear system model. The topology of the various fre- 
quency domain transfer functions in this model is then used to justify the measurement 
and pre-processing steps required to estimate the smoothed scattering impulse response, 
due to a specified double-Gaussian incident plane wave. There are three measurements 
required for each target impulse response estimation. Aside from the actual target 
measurement, a background measurement is needed to subtract out the directlv received 
signal from the transmitting antenna, in addition to that scattered from the various 
chamber surfaces. A third measurement, for a metallic sphere, 1s needed to allow for 
compensation of the transfer function characteristics of the transmitting and receiving 
antennas in the chamber. By effecting a deconvolution comparison with a highly accu- 
rate numerical calculation of the transient scattering response of the sphere target, this 
final measurement also permits the estimation of the target’s response to the same 
excitation as 1s used in the sphere computation. 

Chapter IV considers the details of the actual software that performs the signal ac- 
quisition, signal averaging and deconvolution post-processing of the measured 
waveforms. The acquisition and signal averaging software is a blend of a commercial 
JEEE-488 bus controller (GURU-II) and a custom driver, designed as part of this thesis 
effort. The deconvolution algorithm provides, first of all, graphically displayed iterative 
background subtractions from both the target and sphere measurements to allow com- 
pensation for time drifts. Secondly, the final deconvolution is performed via a variant 


on Riad’s method, as discussed in Chapter III. 
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The validations shown in Chapter V speak for themselves as to the expected meas- 
urement quality obtained using this new transient scattering range. Both low-Q and 
high-Q type targets are considered by the respective comparisons of sphere and thin-wire 
scattering measurements with computed smoothed impulse responses. Further illus- 
trations of the tvpe of measurements to be performed with this range are shown in the 
scale model aircraft backscattering signatures. A simple software modification will allow 
the display of very broadband frequency domain radar cross sections (RCS) for scale 
model targets. 

To quantitatively investigate the fidelity of the measurement system, a noise esti- 
mation of the new implementation was investigated and was modelled. The particular 
signal-to-noise ratio (SNR) of a thick wire scattering measurement was calculated. A 
27 dB SNR of this case demonstrated the relatively high fidelity of the measurement 


system for the case small RCS targets, including silver coated 1/72 scale model aircraft. 


B. FUTURE CONSIDERATIONS 

Even though the demonstrated system performance was excellent, there 1s no such 
thing as too high an SNR in scattering measurements. In this regard, the power- 
bandwidth product has the potential to be enhanced by future efforts. 

The 1 to 7 GHz passband of the scattering system was limited by the transmitting 
GaAS power amplifier. On the other hand, the bandwidth of the HP DPO system can 
go as high as 20 Ghz, while the current antennas have a | to 12 GHz passband. Future 
use of broader bandwidth antennas and power amplifier, having a rated RMS output 
power of 2 watts or more, would increase the overall system performance. More accu- 
rate smoothed impulse responses could be achieved, as well as wider bandwidth RCS 
measurements. 

Because of the limited scope of this thesis, the noise power estimation was onlv an 
approximation. A more comprehensive noise analysis should be carried out in order to 
quantifv and minimize the noise power sources. 

The final consideration is the future investigation of alternate mechanisms for pro- 
viding broadband scattering measurements. One such mature technique is through a 
stepped frequency continuous wave system, employing a network analyzer front-end as 
the coherent receiver, perhaps preceded by a low-noise amplifier to reduce the overall 
system noise figure. A somewhat radical technique, one that 1s currently being investi- 
gated at NPS, is to use a very broadband amplified noise source as the transmitting 


source, while employing the HP DPO as a two-channel correlation receiver. This 


method is akin to a dynamic matched filter, offering enhanced bandwidth, at relatively 
low cost vis-a-vis the stepped frequency setup. Additionally, the potential exists for both 
tactical and strategic employment of such random source systems, either as active 
“pseudo-jammers’ or as totally passive radars relying upon either natural or man-made 


random signal sources for the illumination of targets. 
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APPENDIX A. THE IMPULSE RESPONSE DECONVOLUTION 


The following computer program calculates the Impulse Response from the meas- 
urement data as described in chapter 4. The program is written in RM/FORTRAN. 

The subroutine PLO8 and P11 require the GRAFMATIC routine library. The pro- 
cedure to compile, link and run this source code refers to the RM/FORTRAN User’s 
Guide Version 2.11 (DOS). 


CITT ME Te He FF HH TINE TE HFK FTE KH FTI HAH KEK HNKKKKKKKEE KK 


C THIS IS THE IMPULSE RESPONSE DECONVOLUTION * 
C WRITTEN BY LT. JG SOONPUEN SOMAPEE Feb 88 * 
C ADAPTATION OF ORIGINAL TEKTRONIX PROGRAM 7 
C BY PROF. M.A.MORGAN Circa 1983-85 % 
CS revere dere vodeds Tete deve Tote Ve We ve He Ve He Fe Fe Ve Ve TE TENE Ve Fe Ve Te VE TE TE Te IE FEN Ve Me He He WE TEN Hee 
Cae Hrededededede Teds Ws Fe VeVe ye He FEE Te Ve TENE TeV. FeV TeV HE NC HEHE TER NNER HK KIN 
C WHERE: * 
C DATE IS THE DATE OF MEASUREMENT * 
C DATAB IS THE NUMBER OF DATA BLOCKS % 
C Fl IS THE DUMMY TRING * 
C ID IS THE IDENTIFICATION STRING % 
C K IS THE NUMBER OF POINTS * 
C L9 IS THE RELATIVE SMOOTHING * 
C MAG IS THE MAGNITUDE OF Xo*Ht(f) * 
C NAME4 IS THE NAME OF THE TARGET DATA * 
C N IS THE NUMBER OF POINTS % 
C NUMPT IS THE SAME AS K * 
C NEE IS INDICATEES POWER-OF-2 EXPONENT * 
C WILL BE GENERATE BY PREFFT * 
C SET TO -1 INDICATES ERROR % 
C R IS THE REAL PART OF THE INVERSE * 
C FOURIER TRANSFORM OF Xo*Ht( f£) 7 
C SUBAV IS THE NUMBER OF SUB-AVERAGES * 
C at IS THE TIME INTERVAL * 
C TIME IS THE TIME INTEVAL * 
C Tigi IS THE TITLE NAME FOR PLOTTING % 
C Sey IS THE FREQUENCY DOMAIN OF SSE 
C TYEE IS THE WAVEFORM TYPE 

C W IS THE COMPLEX EXPONENTIAL ARRAY * 
C GENERATED BY PREFFT SUBROUTINE * 
C Xo*Ht(f£) IS THE GAUSSIAN WEIGHTED COSINE ¥ 
GC XS IS THE STRING FOR PLOTTING ROUTINE* 
CG FOR THE HEADER * 
CG DRIVER * 
C Y4 IS TFE CALIBRATION SPHERE * 
C SUBTRACTED WAVEFORM (FREQUENCY ) * 
G * 
C Y5 IS THE TARGET SUBTRACTED WAVEFORM * 
C FREQUENCY DOMAIN ¥e 
C ¥ 
C YSG IS THE COMPUTED SCATTERED ELECTRIC* 
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Y1 
12 
pee 
YS 
YN 
XMI1 
XMA1 


CoieaGa Ga) Ca C2) C2) C2 


FIELD FOR THE CANONICAL TARGET 


Is 
iS 
IS 
IS 
Is 
Is 
IS 


THE 
THE 
THE 
THE 
THE 
THE 
THE 


CALIBRATION MEASUREMENT 
BACKGROUND MEASUREMENT 
TARGET MEASUREMENT 

SAME AS XS& 

YES OR NO STRING 
MINIMUM VALUE 

MAXIMUM VALUE 


CHARACTER*64 TITLE,XS,ID,YS,DATE,TYPE, NAME 
CHeamaciek*64 TYPE ,TYPE2 ,tyPES,F1 
CHARACTER*16 NAME1,NAME2 , NAME3 


CHARACTER*64 ID1,1ID2,I1D3 
CHARACTER*1 


YNS='N' 


YNS,YNS1 
COMPLEX YSC(1025),Y4( 1025) ,Y5( 1025) ,W(1025) 
REAL T,T2 ,TIME 
REAL Y1(1024) ,Y2( 1024) ,Y3( 1024) 


INTEGER N,SUBAV,DATAB , NUMPT 
CARR Hee Te Fe Te Ge Ne Ve Me He Te He Fe Ve TH TeV HVE Fe HE TERK TK KKK IK KICK KICK 


CH eae ede dede Fede Fe He HI VIET He HET TET TEN TCT KE FE TCI EH VE TC HERI I HK 
READ DATA FOR COMPUTED SPHERE THEN 


C 
C 
C 


CALLING SUBROUTINE FFF TO DRAW THE TIME 
SERIES AND SPECTRAL PLOT 


Creve dene dedeve sede de dese de Fede He se ve He de as Te de He Tee Tere Te He Te Tee He Te Fe HE Fe He He Vee ee He Ve He 
l geakeckdesesededescdededevededeaedevededeacsedksedk desesedkeaeskesesdedeskeskesesesdesdeskeaedededeskskesesesesed ' 


WRITE(* ,**) 
WRITE(* ,*) 
WRITE(* ,**) 
WRITE(*,*) 
WRITE(* ,*) 
WRITE(* ,*) 
WRITE(* ,*) 
WRITE(* ,**) 
WRITE(* ,**) 
WRITE(*,*) 
WRITE(*,*) 
WRITE(*,*) 
WRITE(* ,*) 
WRITE(* ,?*) 
WRITE(*,*) 
WRITE(* ,**) 
WRITE(* ,*) 


READ(*,173) NAME 


Toss 
» 


oc 
us 


os. 


t 
t 


Mr 


Vor 
toe 
I oe 
Vx 
Ve 
oe 
Ve 


x 


l Sostevestsleveveveveveveve 


DECONVOLUTION PROGRAM BY DR. M.A. MORGAN 
ADAPTED TO IBM BY LT. JG. SOONPUEN SOMAPEE 


READ 
READ 


OAM E WHE 
> 8) 
td 
> 
Co 


FORMAT OF DATA FILE: 


ID 
TYPE 
DATE 
Pl 


SUB-AVERAGE 
TIME WINDOW 


DATA BLOCK 


( FORMAT 
( FORMAT 
( FORMAT 
( FORMAT 
( FORMAT 
(FORMAT 
( FORMAT 


NUMBER OF POINTS (FORMAT 


Sesetevesevovevetetesededeotedesesvededeteleletedotecdoevciotctclcdekscsdse' 


‘FILE NAME FOR COMPUTED CALIBRATION SPHERE:  ' 


CHARACTER*64 ) 
CHARACTER*64 ) 
CHARACTER*64) 
CHARACTER*64 ) 
INTEGER) 
REAL(F5. 2)) 
INTEGER) 
INTEGER) 


Cr erevededeve dese ve Jeeves veve de vee Fe ve Fe Fe Tee He Ne He He He Te Ve Te VEN He Ve Te He I He Vee He He He Hee He eM He He VEN Hee 


C 
C 
C 
C 
C 
C 
C 
C 


C 


YSC = COMPUTED SPHERE COMPLEX FREQUENCY DOMAIN DATA 


NUMPT = NUMBER OF POINTS IN THE OUTPUT OF THIS ROUTINE 


TIME = OUTPUT TIME INTERVAL 


CALL FFF1( NAME, YSC ,W,NEXPT , TIME) 
Cred rede seve de deve deve de ve He sete Tee He Fee re ve He Te He He eI He He eK Fe NH HFK Te HFK TK HHH Ke KI 


9110 CONTINUE 
Crake dese de se dese dese de deve de de Fede ve ke Fe He ve He eK KI HK CH TICK Fe Fe Te He eK EK NK I KKK 


READ THE DATA FILES FOR CALIBRATE AND BACKGROUND 


47 


Tevekwetcdedededekkksesdcsdededese added wedded kkksdedesdesdekktekesedede keke teste deve tek sede 


x! 
sc! 


*! 


x! 


x! 
ye! 
x! 
a! 
ve! 
*! 
se! 


sc! 


Ca OC) CiO2 G22 C202 C2 C2 co cD 


THEN OVERLAY THESE ONTO THE SAME AXIS. 
PLOT77 ALLOWS THE OPERATOR TO SHIFT ONE WAVEFORM 
PRIOR TO SUBTRACTING THE BACKGROUND. 


THE RESULT FROM SUBTRACTION IS TRANSFORMED TO 
FREQUENCY DOMAIN AND SAVED TO THE COMPLEX ARRAY Y4 


NAME1 = NAME OF THE FILE TO SENT TO PLOT77 ROUTINE 
(CALIBRATE DATA) 

NAME2 = NAME OF THE BACKGROUND DATA 

ID1 = IDENTIFICATION CHRS FOR CALIBRATE 

Y4 = COMPLEX ARRAY OUTPUT (CALIBRATE ) 


Setescacdedcscavcdedeacsdedeseacdesescaksesesdesencakcakeatvacscslescdcdeslkeveaeskcskeatcscsdeseseslcsicsicsedeseseseveveveveve ae vege 


WRITE(* *) l Sgeatedeseckscscskedevedcseveuescddkeavacsdeseaedededeseacveskdedesdededesevesedeseses ' 
’ 


WRITE(* ,*) "FILE NAME FOR CALIBRATION SPHERE: 
READ(*,173) NAME1 

CALL INPPCNAME1,ID1,TYPE1,DATE,F1,SUBAV,T2,DATAB ,NUMPT, Y1) 
WRITE(*,*) l esedeaededcaedesdededeakeaeseacdedeacaeskesesedeededsededededededdovesesedesdesesd | 
WRITE(*,*) "FILE NAME FOR SPHERE BACKGROUND WAVEFORM: ' 
READ(*,173) NAME2 


Crewe dedese devon dese de dete ve Tove te seve te seve Fe ve ve He ve ese ese Ve Fe He Ie We Wee Me Ie ee Wee Fese vere 


CALL INPP(NAME2,ID2,TYPE2,DATE,F1,SUBAV,T2,DATAB,NUMPT, Y2) 
CALL PLO77(Y1,Y2,T2,NUMPT, ID1) 


Cr reaver deeded deveve seve seve deus Fe se Fe ve Tose Teac se se WeNe We We Tere Were Wee He Mee Were we He Ae 
9112 CONTINUE 


WRITE(*,*) V ede desk dese deaedesededeaeaesese dese ak deaeredese desea dese sede sedesedeedede eves | 
3 


WRITE(*,*) ‘FILE NAME FOR TARGET WAVEFORM: 
READ(*,173) NAME3 


Crcdederdededevere deat seve de Tove ve Fe Fe He Te Ne Tete He Ie Ne We Vee We Me Ie TeNe Ve Ve Te He Fe TENE Ne Hee HIE He 


CALL INPP(NAME3, ID3,TYPE3 ,DATE,F1,SUBAV,T2, DATAB ,NUMPT, Y3) 


Comte mt, 


Crcderevervdededevode deve deve sede ye Fede Fe Tete Fede Ne Fe Fe WEF Fe WEI Ne Te VET Te Fe Te Te Ne Fe Ne Fe Ve Fe Me 


WRITE(* ,**) 

WRITE(*,**) "USE PREVIOUS BACKGROUND FOR THE TARGET ? (Y/N): 
READ(**,173) YNS1 

IFC YNS1. EQ. 'Y'. OR. YNS1. EQ. 'y') GO TO 9222 

WRITE(* ,*)  sevesevedededevevedete deve dete dete de OUTER NTR TTI | 
WRITE(*,*) ‘FILE NAME FOR TARGET BACKGROUND WAVEFORM: ' 

READ(* ,173) NAME2 


Creaeaerededevesedede sede de de ae dene He ve devede He vex sexe He dese deve ve dese vere ae ve se seve dese a vee 


CALL INPP(NAME2 , ID2,TYPE2,DATE,F1,SUBAV,T2,DATAB,NUMPT, Y2) 


9222 CONTINUE 


C 
C 
C 
C 
C 
C 
C 
C 


987 


REPEATING SHIFT AND SUBTRACT OPERATIONS ON TARGET AND BACKGROUND 
NAME1= TARGET DATA FILE 

NAME2= BACKGROUND DATA FILE 

OUTPUT 

ID1 = IDENTIFICATION OF THE TARGET 

Y5 = COMPLEX ARRAY OF FREQUENCY DOMAIN TARGET DATA 


Medkscdevededkededededededesedeakcdesdevesedesevesedeskesdesevledke sevese ake sk teste testes de sede de slesese 


CALL PLO77(Y3,Y2,T2,NUMPT,ID3) 

DO 987 I=1,NUMPT 
Y4(I)=CMPLX(Y1(1),0. 0) 
Y5(1)=CMPLX(Y3(1I),0. 0) 

CONTINUE 

MODE=0 
CALL FFT(NUMPT, MODE , TIME ,NEXPT,W,Y4) 
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CALL FFT CNUMPT,MODE , TIME ,NEXPT ,W,Y5) 

WRITE(*,*) OWemeadededevereuecededevere de de seve cesesete sete sete eh IKI | 

CALL DECON3(NUMPT,YSC,Y4,Y5,ID3,TYPE3 ,DATE,F1,SUBAV,T2,DATAB) 
WRITE(* ,*) UA ET AeA Ar ok Ar Gideon CR hc oe pe ate ae aie aR ek i get kek a aie id el 
WRITE(*,*) ‘WOULD YOU LIKE TO PROCESS ANOTHER TARGET ? (Y/N): ' 
ReaD ( 173) YNS 


IF(YNS.EQ. 'Y' .OR. YNS.EQ.'y' ) GO TO 9112 
WRITE(* % ) " Sevedevedevedevededevevededetededetededesdedeleskdedledetedek ' 


b 
WRITE(*,*) '* eeeee eeeeece e *! 
WRITE(*,*) '* eee eee eeee ve! 
WRITE(*,*) '‘eeeeeeeeeeeeeeeeeeeeeeeeeeeee *' 
WRITE(*,*) ' eeeececeeeeeceececeeeeeeeee iY 
WRITE(*,*) ' eeeeeeeceeeeeeeeeeeeeee ve! 
WRITE(*,*) " ededededededevedetdedetededetetenedletekietesecntedetictedete | 


’ ' ’ teat U U ' 
WRITE(*,*) TeteteFeFeHs Fete de Ke Tove FeV TET FeV Te Te He PICNIC NN RK 


173. + FORMAT(A) 
»lOP 
END 


SUBROUTINE INPP( NAME, ID,TYPE,DATE,F1,SUBAV,T2,DATAB,NUMPT, Y) 


Creede dese Te seco ee TORT TERETE TS TOT E TET TS BETS TOTES ETE TERE DETER TELE TR EVER ETE EVER ARE TERE ERE RRER ETRE 


C 

é INPUT NAME 

C OUTPUT ID,TYPE,DATE,F1,SUBAV,T2,DATAB,NUMPT, Y 

C NAME = CHRS NAME OF THE DATA FILE 

C ID = CHRS TARGET IDENTIFICATION 

C TYPE = CHRS WAVE FORM IDENTIFICATION 

C DATE = CHRS DATE 

C INTEGER 

C SUBAV = NUMBER OF SUB-AVERAGES USED TO ACQUIRE THE DATA 
C DATAB = DATA BLOCK 

C Fl = DUMMY 

C NUMPT = NUMBER OF POINTS 

C REAL 

C 2 = TIME WINDOW 

C ag = Y ARRAY OUTPUT 

Crerevedesevedete See Tete lene ce te Fe OTe Te ae TOTS TETAS TE TEDET ETE LODEDE LEDC TE PED EPEDE TERED PERE TERE TERE TRE T ETON BEER EG 


CHARACTER*16 NAME 
CHARACTER*64 ID,F1 
CHARACTER*64 TYPE ,DATE 
REAL*4 Y¥(1500) 

INTEGER SUBAV,DATAB,NUMPT,K 
PI=3. 141592654 
Kz} 
OPEN (UNIT=10,FILE=NAME , STATUS='OLD' ) 
© READ THE HEADER 
READ(10,5) ID 
READ( 10,5) TYPE 
READ(10,5) DATE 
READ(10,5) Fl 
READ( 10,2030) SUBAV 
READ( 10,2010) T2 
READ( 10,2030) DATAB 
READ( 10,2030) NUMPT 
10 READ(10,2000,END=20) Y(K) 


49 





K=K+1 
1? GOTO 10 
20 CLOSE( 10) 
2000 FORMATC E12. 6) 
2010 FORMATCFS. 2) 
2030 FORMAT( I5) 
2, FORMAT(A) 
K=Khea 
IFCK. EQ. NUMPT Seamer sa00 
WRITE(*,*) ‘CHECK THE DATA - # POINTS MISMATCH' 
STOP 
3000 RETURN 
END 


SUBROUTINE FFF1C NAME,YY,W,NEXPT,TIME) 
Crededededededevedecevede se dete de deve Fee dete ve seve Tete Hee Tete Le FoF TEN HIKER RE 


C WHERE * 
C DATE IS THE DATE OF MEASUREMENT * 
C DATAB IS THE NUMBER OF DATA BLOCKS * 
C Fl IS THE DUMMY STRING * 
C ID IS THE IDENTIFICATION STRING * 
C K IS THE NUMBER OF POINT * 
C L9 IS THE RELATIVE SMOOTHING * 
C MODE IS EQUAL TO 0 MEAN FORWARD * 
C FOURIER TRANSFORM 

C MAG IS THE MAGNITUDE OF Xo*Ht(f) * 
C NAME4 IS THE NAME OF THE TARGET DATA * 
C N IS THE NUMBER OF POINTS * 
C NUMPT IS THE SAME AS K * 
C NEXPT IS THE POWER-OF-2 EXPONENT %* 
C WILL BE GENERATED BY PREFFT * 
C SET LO. =i). INDICATES ERROR * 
C R IS THE REAL PART OF THE INVERSE * 
C FOURIER TRANSFORM OF Xo*Ht(f) * 
C SUBAV IS THE NUMBER OF SUB AVERAGE * 
C i IS THE TIME INTERVAL * 
C TIME IS THE TIME INTEVAL * 
C shine IS THE TITLE NAME FOR PLOTTING ¥ 
C ALVA IS THE FREQUENCY DOMAIN OF Xo*H( f)* 
C dhsélele IS THE WAVEFORM TYPE * 
C W IS THE COMPLEX EXPONENTIAL ARRAY * 
C GENERATE BY PREFFT SUBROUTINE * 
C Xo*Ht(f) IS THE GAUSSIAN WEIGHTED COSINE % 
C XS IS THE STRING FOR PLOTTING ROUTINE* 
C FOR THE HEADER * 
C DRIVER * 
C Y4 IS THE CALIBRATION SPHERE * 
C SUBTRACTED WAVEFORM (FREQUENCY) * 
Gc nS 
C x2 IS THE TARGET SUBTRACTED WAVEFORM * 
C FREQUENCY DOMAIN * 
G ve 
C YSC IS THE COMPUTED SCATTERED ELECTRIC* 
C FIELD FOR THE CANONICAL TARGET * 
C yell IS THE CALIBRATION MEASUREMENT % 
C TZ IS THE BACKGROUND MEASUREMENT ¥ 
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= 


C NES IS THE TARGET MEASURED % 
C Yo IS THE SAME AS XS * 
C YN to THEY YBS OR@NO STRING * 
C MEEK * 
C I IEA IS THE MINIMUM VALUE * 
C XMA1 IS THE MAXIMUM VALUE * 
Ctedede deve deve de deve He Fe Fe Weve Fede He TFG He Hee Te TE HET ITE KE TEN TE IT RE RET 


CHARACTER*16 NAME 
CHARACTER*64 ID,TYPE,DATE,TITLE,XS,YS,F1 
COMPLEX W(1025),YY( 1025) 
REAL T 
REAL Y( 1024) ,TIME,T2,XMINS ,MAG( 1024) , XMAG 
INTEGER NUMPT, DATAB ,SUBAV,MODE,N 
Credo dese ded de devede te Fete Fee F8 Fe Fe eT HE He TCT Fe HET KE TET TET TN ETEK RK IK 
CALL INPP( NAME,ID,TYPE,DATE,F1,SUBAV,T2,DATAB,NUMPT,Y) 
TIME=T2/(FLOAT( NUMPT-1) )*1E-9 
XMINS=0. 0 
=NUMPT 
CALL PRETIM(TITLE,XS,YS) 
CALL P11(TITLE,N,XMINS,T2,Y,XS,ID,YS) 
C SET MODE EQUA: TO 0 TO FIND FOURIER SERIES 
MODE=0 
CALL PREFFT (N,MODE,NEXPT,W) 
DO 1234 I =1 ,NUMPT 
YY(I)=CMPLX( Y(1I),0.0 ) 


1234 CONTINUE 
T=TIME 
CALL FFT (N,MODE,T,NEXPT,W,YY) 

C CHANGING THE OUTPUT YY TO MAGNITUDE AND PHASE 
N=NUMPT 


CALL PREFRE(TITLE,XS, YS) 

CALL MAGNITUD( YY,N, MAG) 

XMAG=1. 0/TIME/2. O* 1E-9 

CALL P11(TITLE,N,XMINS,XMAG,MAG, XS, ID, YS) 


RETURN 

END 
Cre dedeveve tess teva revered de Ted de TeRe TTT Bee TEE DETTE DOVE TE TCT OLE TE TO VE DESERTS TERETE AIR EET 
@ * 
C FFT PROGRAM FROM THE DIGITAL SPECTRAL ANALYSIS 7 
C WITH APLICATIONS PRENTICE-HALL, INC.. * 
C ENGLEWOOD CLIFF,NEW JERSEY 07632 1987 7 
c S. LAWRENCE MARPLE * 
c J 
Cred dedkde deve rede de Fede sevede de FX He TON TT KK Fe eFC I TR RT TCR RCC RR RRR TK RET 

SUBROUTINE PREFFT (N,MODE ,NEXP,W) 

Input Parameters: 
N - Number of data samples to be processed (integer-must be a 


power of two) 
MODE - Set to 0 for discrete-time Fourier series (Eq. 2.C.1) or 
1 for inverse (Eq. 2.C. 2) 


Output Parameters: 


G39Ga.G) C2 CaeG2 CC) C) Ga) 


5] 


C) GC) Cl Gases) G2 C2 C2 


10 


Lp 


20 


NEXP - Indicates power-of-2 exponent such that N=2**NEXP . 
Will be set to -1 to indicate error condition if N 
is not a power of 2 (this integer used by sub. FFT) 

W - Complex exponential array 


Notes: 
External array W must be dimensioned .GE. N by calling program. 


COMPLEX W( 1024) ,C1,C2 
NEXP=1 
NT=2%**NEXP 
IF (NT .GE. N) GO TO 10 
NEXP=NEXP+1 
GO TO 5 
IF (NT .EQ. N) GOTO 15 
NEXP=-1 
RETURN 
S=8. *ATAN(1. ) /FLOAT(NT) 
C1=CMPLX(COS(S),-SIN(S)) 
IF (MODE .NE. 0) C1=CONJG(C1) 
C2=(1. ,0.) 
DO 20 K=1,NT 

W(K)=C2 

C2=C2*Cl 
RETURN 
END 


Crcdere dese ese deve deve de Hee ve Te eye FoF Fe MEIC TENE Te Me HE HE VEE ETE NE TE Fe ETE He Te FEE TE I ee EE 


GleCeaGa Cu Gn cuca Cl) Cy C2 C2 C2 CIC? CC) © 


SUBROUTINE FFT (N,MODE ,T,NEXP,W,X) 
Input Parameters: 


N,MODE,NEXP,W - See parameter list for subroutine PREFFT 
T - Sample interval in seconds 
X - Array of N complex data samples, X(1) to X(N) 


Output Parameters: 


X - N complex transform values replace original data samples 
indexed from k=1 to k=N, representing the frequencies 
(k-1)/NT hertz 


Notes: 
External array X must be dimensioned .GE. N by calling program. 


COMPLEX X( 1024) ,W(1024),C1,C2 
MM=1 
LL=N 
DO 70 K=1,NEXP 

NN=LL/2 

JJ=MM+1 

DO 40 I=1,N,LL 

KK=I-+NN 
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C1=X( 1 )+X(CKK) 
X(KK)=X( 1) -X(CKK) 
40 X(I)=C1 
IF (NN .EQ. 1) GO TO 70 
DO 60 J=2,NN 
C2=W( JJ) 
DO 50 I=J,N,LL 
KK=I+NN 
C1=X(1)+X(KK) 
X(KK)=( XC 1) -X(KK) )*C2 


. 50 X(-c1 
60 JJ=JI+MM 
LL=NN 
MM=MM** 2 
70 CONTINUE 
NV2=N/2 
NM1=N~1 
J=1 


DO 90 I=1,NM1 
PeolieaGk. J) GO TO 80 


C1=X(J) 
X(J)=X(1) 
X(I)=C1 
80 K=NV2 
85 iEkeNGE. J) GO TO 90 
=J-K 
K=K/2 
GO TO 85 
90 J=J+K 
IF (MODE .EQ. 0) S=T 
IF (MODE .NE. 0) S=1. /(T*FLOAT(N)) 
DO 100 I=1,N 
100 X(I)=X(1)*S 
RETURN 
END 


Cr aeH Tete AK PR TT eR IR TKR RK RRR ER ERE RI 


C 
SUBROUTINE MAGNITUD (Y,NUMPT, MAG) 


Cr dere tetedede deaede Feeds Fete HR TRI RTE RINK IK RRR TERRE KERR 
COMPLEX Y( 1024) ,H 
REAL*4& MAG( i024) 
INTEGER N,K,NP,NUMPT 
DO 10 I=1 ,NUMPT/2 
H=Y( I) 
MAG( I )=CABS(H) 
10 CONTINUE 
NUMPT=NUMPT / 2 
RETURN 
END 


Crete dedede vedas dedede Heese dese Fe dede de KK RICE RC RE RTE RI 


SUBROUTINE PRETIM(TITLE,XS,YS) 
CHARACTER*64 TITLE,XS,YS 


a0 


TITLE=' TIME SERIES PLOT' 


XS  ='X-AXIS TIME IN NS ' 

YS  ='Y-AXIS MAGNITUDE IN VOLT' 
RETURN 

END 


Cre Hae Tete de ve Te Tete te Fo7e He ve Fe He Fe Te Fe se He Fee Te Fe Te I Te Fe Ie Fe Te Fe Ne Fee He He He Fe We Ne He WH 


SUBROUTINE PREFRE( TITLE, XS,YS) 
CHARACTER*64 TITLE,XS,YS 
TITLE='SPECTRAL PLOT' 
XS='X-AXIS FREQUENCY IN GHZ' 
YS='Y-AXIS MAGNITUDE IN VOLT' 
RETURN 

END. 


SUBROUTINE IFFF(N,T,X) 
Creve eve Fese as ye deve Fedewe He Fe Fee Ve Fe Te Ve HE Ie Me He Te WIE NE Fe ve Fe Te WEN Fe IE HEE Te Te NE Te TEN TeV eNe We Te Fe Fe He Fe Fe Fe Fe Fe 
COMPLEX W(1025) ,X( 1025) 
REAL T 
INTEGER MODE 
MODE=1 
CALL PREFFT(N, MODE ,NEXP,W) 
CALL FFTCN,MODE,T,NEXP,W,X) 
RETURN 
END 


Cr Wasse dew as se de Weve re seve se de Ye Weg ae Ve Fe Ve He ae Fee Ve He Ve VEN He Fe Ve IE NE IE Te MEW IE Te He WEE Ve Fe Ie Me Wee Fe Fe Ie Te Ve Ie Fe Ie He Ie He Ie WN Tee 


SUBROUTINE DECON3(K,YSC,Y4,Y5,ID,TYPE,DATE,F1,SUBAV,T2,DATAB)} 


¢ 


Creve detec ce eM rs sede deve te He se Te Tove Be VET Tose Fe Te Ts He Te TENE Fe Fe Te Ve Fe Tee Fe Te Fe VET VEE TEN EN HF 


& COMPUTING THE TARGET RESPONSE TO THE INCIDENT * 
& WAVESHAPE FOR THE COMPUTED SPHERE. USING * 
C RIAD'S METHOD OF OPTIMAL DECONVOLUTION. * 
C aC 
& Xo*Ht( f)=Y5*CONJCY4)*Ysc/(Y4*CONJ(Y4)+L) * 
CG a 
C WHERE %* 
c Xo(f) TS THE INCIDENT WAVEFORM * 
'@ ¥ 
C Ht(£) IS THE TARGET TRANSFER FUNCTION % 
C es 
C Y4 IS THE CALIBRATION SPHERE ¥ 
C SUBTRACTED WAVEFORM (FREQUENCY ) * 
@ * 
C y5 IS THE TARGET SUBTRACTED WAVEFORM * 
C FREQUENCY DOMAIN * 
C ve 
G Yse IS THE GOMPUTED SCATTERED ELECTRIC ¥ 
CG FIELD FOR THE CANONICAL TARGET ¥ 
C al IS THE CALIBRATION MEASUREMENT * 
C Y2 IS THE BACKGROUND MEASUREMENT * 
C nes IS TrE TARGET MEASURED - * 
C K IS THE NUMBER OF POINT ve 
G NUMPT IS THE SAME AS K 7% 
C R IS THE REAL PART OF THE INVERSE ¥ 
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C FOURIER TRANSFORM OF Xo*Ht(f) * 
C MAG IS THE MAGNITUDE OF Xo*Ht(f) * 
C ID IS THE IDENTIFICATION STRING * 
C SUBAV IS THE NUMBER OF SUB AVERAGE * 
C DATAB IS THE NUMBER OF DATA BLOCK * 
C Fl IS THE DUMMY STRING * 
C NAME4 IS THE NAME OF THE TARGET DATA * 
C XS IS THE STRING FOR PLOTTING ROUTINE * 
C FOR THE HEADER * 
C YS IS THE SAME AS XS * 
C TITLE IS THE TITLE NAME FOR PLOTTING ¥ 
C YN IS THE YES OR NO STRING * 
C L9 IS THE SMOOTHING PARAMETER * 
C Jee IS THE FREQUENCY DOMAIN OF Xo*H(f) * 
C ie E IS THE WAVEFORM TYPE * 
C XMI1 IS THE MINIMUM VALUE * 
C XMA1 IS THE MAXIMUM VALUE * 
CG DATE IS THE DATE OF MEASUREMENT * 
Cw Hae sede Neve Fe ve ve se We Fe ae ae Wee Te Wee Hee Ve Fe Ne Fe Hee Fee Ie Fee Fee Ve Hee WME He He Ve He Fee He Fe He Hee 

COMPLEX V5 C1025e r4( 1025) -YSC( 1024) ,TUTC 1024) ,SUM 

COMPLEX | OE 29 I 

REAL R( 1024) ,L9 ,XMA1,XMI1,MAG( 1024) 

REAL I ts 9 Ape heed hae 

INTEGER K ,NUMPT , SUBAV , DATAB 


CHARACTER*16 NAME4 
CHARACTER*64 ID,YS,XS,TITLE,TYPE,DATE,F1 
CHARACTER*1 YN 
Crevededede deck deve stese de ss deve Fe se ve He ve ve tose ae He se Teak Je Ne He ae Hee Te He Heese Ne ve We eae ve He He 
T=T2 /FLOAT(K-1)*1E-9 
SUM=(0. 0,0. 0) 
DO 111 I=1,K 
SUM=SUM+Y4(1)*CONJG(Y4(I)) 
111 CONTINUE 
438 WRITE(*,*)  Sevevedededededededesedetededetesesesetede eK TRIER IERIE | 
WRITE(*,*) ' ENTER RELATIVE SMOOTHING PARAMETER : 
READ (*,*) L9 
L1=CMPLX(L9 ,0. 0) 
L2=SUM/CMPLX( FLOAT(K+1) ,0. 0)*L1 
DO 1 I=1,K 
TUT( I )=Y5(1)*CONJG(Y4(1))*YSC(1)/CY4(1)*CONJG(Y4( 1) )+L2) 
i CONTINUE 
TITLE='SPECTRAL PLOT' 
YS=' DECONVOLUTION' 
XS='X-AXIS FREQUENCY IN GHZ' 
NUMPT=K 
C everest sede Te Tee Te eH TTS TTS TE CTE PTE ET TERETE OLE TER TELLER RRR RE RERERREK 


CALL MAGNITUD( TUT , NUMPT , MAG) 

Ce detede secede dele Feed Fede Fede TOG TEHe TE TS CTE TE TRIES TTT TET RIK TREE ERRECER 
XMI1=0. 0 
XMA1=1. 0/T/2. O*1E-9 

Credeacdedevdedede sealed deve dere Jevese Heese eve ae He He Fee He Fete Fe He Te ete He eK He Hee HK Te Fee Ke 


CALL P12(TITLE,NUMPT,XMI1,XMA1,MAG,XS,ID,YS,L9) 
Crerevedededevede teak ese He ve de se ve ve Fe ae ae ae ve Ve Fee We Ne Fe Te Fe He Fee Ie He Fe Ne Fe ae He Fe He Hee He He eK 


GALL IFFF(K;T,TUT) 
Crvavaededesedeaedeaese ve devede ve ae ve ved ve ae ve ae ve se He ae He see Ve se ae Tes He He ve Te Fee He Ne te He He Fe HK 


DO 222 I=1,K 
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Cee 


RC I)=REALC TUT( I) ) 
CONTINUE 


Creve dere aede ve de dese evs HF He Ie He Tee WE Te We IE Te ae We Me Ve IE Te se IE Me Fe Ie Vee Ie Ve IE Tee We Te Ie Ie ae He Ves 


XMI1=0. 0 

XMA1=T2 

XS='X-AXIS TIME IN NSEC' 
TITLE='TIME SERIES PLOT' 
YS='DECONVOLUTION' 


Crevede de dede de dede ve We Nese we Fe We 7s Tes He Te Beas Te He He Fe He IE Ie He ee He Me He He Tee Hee He Fe He Ie Fe Feve He 


CALL P12(TITLE,K,XMI1,XMA1,R,XS,ID,YS,L9) 


Crededededededede de de deve dened ve ded desea ss deve Heese se Fee Fe Neve dese Tete Wea Te se Nee Ve Hee He 


YN='N' 

WRITE(* ,**)  Jeavedkdeveskdeaededededeaeakaeakaeded dedede de ded ded dedede dete ded dese Head eed de desea dese sede | 
WRITE(* ,**) ‘WOULD YOU LIKE TO SAVE THE DECONVOLUTION RESULT ?: 
READ(* ,999) YN 

IF ( YN. EQ. Y .OR. YN. EQ. y )°GOmnge eo 

GOTO 4388 


Cree dede dete ds ve Hue He ve He We ss Vege Fe Te We He HEN eT Tose Wee He Te He Tee Te Te ere seve Te Ke weNere Ww Here 


4376 


4388 


992 


wT ote etectealon? Pelee! ’ ’ ’ ve? ’ t 
WR PPE se ; % ) sevcocdcveaveveosvenvovvedke sea dkeocakvocseskvskeakakakovovakakeakscake kak kako 


WRITE(**,*)' ENTER FILE NAME FOR THE TIME DOMAIN DATA:  ' 
READ(*,999) NAMEG 

CALL OUT(NAME4 , ID, TYPE ,DATE,F1,SUBAV,T2,DATAB,K,R) 

YN='Y 

WRITE (ceo) l eaeaeaeaeseaedeskdedeasaedeaesk deakeaeaeeseaedededededede ae deaedevedkaedededeaeaededkeakeded desde sen | 
WRITE (*,%*) "REPEAT WITH NEW SMOOTHING PARAMETER ? (Y/N): ' 
READC* ,999) YN 

IF ( YN. EQ. 'Y'.OR. YN. EQ. 'y') GO TO 438 

FORMAT( A) 

RETURN 

END 


C were aeaede dese ce He dev ae ae Fs He Tone dese eae Teds ve te te ve sede ae te He dee Fe vee Fe ye He eae Tee Ie He He He 


PE Se Sp al Se 


eee) Cela cunG2C)) Ca'Ga Co CaG7 C2 


F = ARRAY OF Y VARIABLE 

XS= $ X-SCALE (STRING) 

YS= $ Y-SCALE (STRING) 

ID= IDENTIFICATION OF THIS GRAPH 


SUBROUTINE P1li(TITLE,NPTS, XMINN, XMAXX ,FFF,XS,ID,YS) 


RM-FORTRAN Subroutine to Plot a Solid Line with TITLE 
Based on Program PLOT. FOR. Dec 1987 by M.A. Morgan 


INPUT DATA FORMAT 


TITLE - 64 Space Title 
NPTS - # Data Points 
XMIN - Real Min X value 
XMAX - Real Max X value 
F(N) - Input Data Array 
THC - Hardcopy (1=Yes) 


CHARACTER*1 YN,DUM 
CHARACTER*64 TITLE,XS,ID,YS 
REAL X(1025),F( 1025) ,XMIN, XMAX, FFF( 1025) , XMINN, XMAXX 
INTEGER*2 N,JROW,JCOL, ISYMBL, ITYPE,CYAN,WHITE, YELLOW 
LHC=0 
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WRITE(*,*) ' WOULD YOU LIKE TO HAVE A HARD COPY: ' 
READ(*,100) YN 
IF(YN. EQ. 'Y'. OR. YN. EQ. 'y') IHC=1 


WHITE=7 
CYAN=11 
DO 1111 J=1,NPTS 
FC I)=FFFC(T) 
ic lbaia CONTINUE 
XMIN=XMINN 
XMAX=XMAXX 
YELLOW=14 
N=NPTS 
DX=( XMAX-XMIN)/(NPTS-1. 0) 
FMIN=0. 0 
FMAX=0. 0 


DO 22 K=1,NPTS 
X(K)=XMIN+(K-1.0)*DX 
IF(F(K). LT. FMIN) FMIN=F(K) 
IFC F(K).GT. FMAX) FMAX=F(K) 
2 CONTINUE 
IF( FMIN. GT. 0.0) FMIN=0. 0 
IF( FMAX. LT. 0.0) FMAX=0. 0 
Computing Scale Factors for Vertical Axis 
ABSMIN=ABS(FMIN) 
ABSMAX=ABS( FMAX) 
YMAX=AMAX1( ABSMIN, ABSMAX) 
NSCL=INT( LOG10( YMAX) ) 

IF (YMAX. LT. 1.0) NSCL=NSCL-1 
YSCL=10. **NSCL 
FMIN=FMIN/YSCL 
FMAX=FMAX/YSCL 

ABSMIN=ABSMIN/YSCL 
ABSMAX=ABSMAX/YSCL 
DO 33 K=1,NPTS 


oe FC K)=FC(K)/YSCL 
YMIN=0. 0 
TECEMINGEO=0.0;) GO 10.37 
25 YMIN=YMIN+0. 5 


IFC ABSMIN. GT. YMIN) GO TO 35 
YMIN=YMIN*FMIN/ABSMIN 


37 CONTINUE 
YMAX=0. 0 
IF(FMAX. EQ. 0.0) GO TO 41 
39 YMAX=YMAX+0. 5 


IFC ABSMAX. GT. YMAX) GO TO 39 
YMAX=YMAX** FMAX / ABSMAX 
41 CONTINUE 
CALL QBEEP 
CALL QSMODE( 2) 
WRITE(*,*) ‘Press RET to View Plot...RET Again to Clear Screen' 
IFC IHC.NE.1) GO TO 42 
WRITE(*,*) ‘A Hardcopy Will Be Made...Check that Printer is ON' 
42 CONTINUE 
READ(*,100) DUM 
C Calling GRAFMATIC Routines and Plotting Solid Line Graph 
ITYPE=1 


2) 


ISYMBL=-2 
CALL QSMODE( 16) 
CALL QPTXT(64,TITLE, YELLOW, 16,23) 
CALL QPTXT(64,ID,CYAN,16,19) 
CALL QPTXT(64,XS,CYAN, 16,20) 
CALL QPTXT(64,YS,CYAN, 16,21) 
CALL QPLOT(110,540,20, 240, XMIN, XMAX, YMIN, YMAX,0. ,0. ,0,1. ,1. 5) 
CALL QSETUP(0,CYAN, ISYMBL,CYAN) 
XMAJOR=XMAX/5. 0 
CALL QXAXIS( XMIN, XMAX, XMAJOR,1,1,2) 
YMAJOR=1. 0 
CALL QYAXIS( YMIN, YMAX, YMAJOR,1,1,2) 
JROW=14+( ABS( YMIN) / (ABS( YMAX)+ABS(YMIN) ) )*220 
JCOL=80 -430*XMIN/( XMAX-XMIN) 
CALL QGTXT(3,'0.0' ,WHITE,JCOL, JROW,0) 
GALL, ORDXTClees , YEINOW, 5. 16) 
CALL OPRTXI.G1., 'c' VLOW 5, 15) 
NOL QPTXT(1,'a' ,YELLOW,5, 14) 
CALL QPTXT(1,'1' , YELLOW ,5,13) 
CALL QPTXT(1,'e',YELLOW,5,12) 
GALL OPTXT(1, X’ , YEUMOW 5.10) 
CALL QPTXT(2,'10' , YELLOW,4,8) 
CALL QPTXT(1,'*', YELLOW, 5,7) 
CALL QPTXTG1, *"FYREG Wes 6) 
CALL QCMOV(0,6) 
WRITE(*,150) NSCL 
CALL QTABL(ITYPE,N,X,F) 
READ(*,100) DUM 
IF(IHC.NE.1) GO TO 44 
CALL PRTSC 
44 CONTINUE 
CALL QSMODE(2) 
100 FORMAT(A) 
150 FORMAT( 4X ,13) 
RETURN 
END 


Cis dedete ted Tere Fete cede Teds TOT PETE TET TESS TOE SOTCTS TOT LOSE TODOS TE TEESE TORE TER ETO TCE 
THIS SUBROUTINE PLOTS TWO GRAPHS ON THE SAME AXIS 


F 
Fl 


C 

C 

C 

C Y1 ARRAY VARIABLE ( BLUE) 

C Y2 ARRAY VARIABLE (GREEN) 

C AFTER PLOT THIS PROGRAM BE AWARE OF THE F1,F 
C IT WILL SENT BACK A SCALE F,F1 
C 

C 

C 

C 

C 


XS = $ X-SCALE 

YS = $ Y-SCALE 

ID IDENTIFICATION 
SevekKkKevekckiekeenskiKhKkKeekidthksKkekkkkRRRSichevek 


SUBROUTINE PLO8( TITLE ,NPTS , XMINN, XMAXX, PUM, XS,YS,PUM1, ID) 


RM-FORTRAN Subroutine to Plot a Solid Line with TITLE 
Based on Program PLOT. FOR. Dec 1987 by M.A. Morgan 


CjaconG. CiuG3 


INPUT DATA FORMAT 
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TITLE - 64 Space Title 
NPTS - # Data Points 
XMIN - Real Min X value 
XMAX - Real Max X value 
FCN) - Input Data Array 
THC - Hardcopy (1=Yes) 


Geen Ga G2 C2°C2 C2 C2 


CHARACTER*1 YN,DUM 
CHARACTER*64 TITLE,XS,ID,YS 
REAL X( 1025) ,F( 1025) ,F iC 1025) ,PUM( 1025) ,PUM1( 1025) , XMINN, XMAXX 
INTEGER*2 N,JROW,JCOL,ISYMBL,ITYPE , CYAN ,WHITE , YELLOW , GREEN 
Created dese de sede dese de seve Fev Teese Feve Fe Fe ve Fe se se Fe se Ne He ese He se Fe Ne Ve Tee Te HEN Tee ve Tee He 

WRITE(*,*) ‘WOULD YOU LIKE TO HAVE A HARD COPY ?: ' 
READ(*,100) YN 
IHC=0 
ieG@vilenoe Y .OR. YN-EOQ. y') IHC=1 
DO 159 I=1,NPTS 

FC I)=PUM( TI) 

F1( I )=PUM1(1) 

159 CONTINUE 


XMIN=XMINN 
XMAX=XMAXX 
Creaeavdeaeakaeaecededeak ae seve aeveve.. dove de deae ve dese ae ae ve se ae seve se ae ae ae ae se ae ae ae ae Me se ae 
WHITE=7 
CYAN=11 
YELLOW=14 
GREEN=2 
N=NPTS 
DX=( XMAX-XMIN)/(NPTS~-1. 0) 
FMIN=0. 0 
FMAX=0. 0 
FMIN1=0. 0 
FMAX1=0. 0 


DO 22 K=1,NPTS 
X(K)=XMIN+(K-1. 0)*DX 
IF(FCK). LT. FMIN) FMIN=F(K) 
IF(FCK).GT. FMAX) FMAX=F(K) 
IF(F1(K). LT. FMIN1) FMIN1=F1(K) 
IF(F1(K). GT. FMAX1) FMAX1=F1(K) 
ZZ CONTINUE 
FMIN=AMIN1I( FMIN,FMIN1) 
FMAX=AMAX1( FMAX , FMAX1) 
IFC FMIN. GT. 0.0) FMIN=0. 0 
IF( FMAX. LT. 0.0) FMAX=0. 0 
C Computing Scale Factors for Vertical Axis 
ABSMIN=ABS(FMIN) 
ABSMAX=ABS( FMAX) 
YMAX=AMAX1( ABSMIN , ABSMAX) 
NSCL=INT( LOG10( YMAX) ) 
IF (YMAX. LT. 1.0) NSCL=NSCL-1 
YSCL=10. **NSCL 
FMIN=FMIN/YSCL 
FMAX=FMAX/YSCL 
ABSMIN=ABSMIN/YSCL 
ABSMAX=ABSMAX/YSCL 
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DO 33 K=1,NPTS 
FIC K)=F1(K)/YSCL 


3) FC K)=FCK)/YSCL 
YMIN=0. 0 
IF(FMIN. EQ. 0.0) GO TO 37 
pe. YMIN=YMIN+0. 5 


IFC ABSMIN. GT. YMIN) GO TO 35 
YMIN=YMIN*FMIN/ABSMIN 


7 CONTINUE 
YMAX=0. 0 
IFC FMAX. EQ.0.0) GO TO 41 
39 YMAX=YMAA+0. 5 


IFC ABSMAX. GT. YMAX) GO TO 39 
YMAX=YMAX*FMAX/ ABSMAX 
41 CONTINUE 
CALL QBEEP 
CALL QSMODE(2) 
WRITE(*,*) ‘Press RET to View Plot...RET Again to Clear Screen' 
IFCIHC. NE. 1) GO TO 42 
WRITE(*,*) ‘A Hardcopy Will Be Made...Check that Printer is ON' 
42 CONTINUE 
READ(*,100) DUM 
C Calling GRAFMATIC Routines and Plotting Solid Line Graph 
ITYPE=1 
ISYMBL=-2 
CALL QSMODE( 16) 
CALL QPTXT(64, TITLE, YELLOW, 16,23) 
CALL QPTXT(64,ID,GREEN, 16,19) 
CALL QPTXT(64,XS ,WHITE,16,20) 
CALL QPTXT(64,YS,CYAN,16,21) 
CALL QPLOT(110,540,20,240,XMIN, XMAX , YMIN, YMAX,0. ,O. , 0) eee 
CALL QSETUP(0,CYAN, ISYMBL,CYAN) 
XMAJOR=XMAX/5. 0 
CALL QXAXIS( XMIN, XMAX ,XMAJOR,1,1,2) 
a YMAJOR=1. 0 
CALL QYAXISCYMIN, YMAX, YMAJOR,1,1,2) 
JROW=14+( ABSC YMIN) /CABS( YMAX)+ABSC YMIN) ) )*220 
JCOL=80-430**XMIN/ CXMAX-XMIN) 
CALL QGTXT(3, 0.0) {WHITE IGOR, JROWAON 
CALL OQPTXT@],'s ,YEmnOW,5, 16) 
CALL QPTXT(1, c. | YELLOW 5.15) 
CALL QPTXT(1, a JYELLEOW, 5, 14) 
CALL OPTXTG1,.0 1 ye Ee@n >, 13) 
GALE L@PIhGl, @ YEON e5. l2) 
CALL QPTXT(1,'X‘ , YELLOW,5,10) 
GALL OPIXT 2h 10 we LOW 4,8) 
GALL GPIXT(1, * ,YELEOW,5,7) 
CALL QPTXT(1,'*' , YELLOW,5,6) 
CALL QCMOV(0,6) 
WRITE(*,150) NSCL 
CALL QTABLCITYPE ,N,X,F) 
Crevedede deve dete Fede secede ve eR TER TTL TTR TEE TETRIS 


G PLOT ANOTHER GRAPH 
Cre tete tele Feve cede te ede STs TENET TE TOTES EERE BEARER ERE 
GREEN=10 


ISYMBL=-2 
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ITYPE=1 
CALL QSETUP(0,GREEN, ISYMBL, GREEN) 
CALL QTABL(ITYPE,N,X,F1) 
READ(**, 100) DUM 
IF(IHC.NE.1) GO TO 44 
CALL PRTSC 
44 CONTINUE 
CALL QSMODE(2) 
100 FORMAT( A) 
150 FORMAT( 4X, 13) 
RETURN 
END 


Cte desede sedeue ve ve sede sede se seve ve He Fe Te deve Hee Fe Te vese Te Fe Tose VEN Fe Fe Fe He Fe Fe Fee He Fe We Fee HK 


C F = ARRAY OF Y VARIABLE 
C XS= $ X-SCALE (STRING) 
C YS= $ Y-SCALE (STRING) 
C ID= IDENTIFICATION OF THIS GRAPH 
C 

SUBROUTINE P12( TITLE,NPTS, XMINN, XMAXX,FFF,XS,ID,YS,L9) 
a 
C RM-FORTRAN Subroutine to Plot a Solid Line with TITLE 
C ‘Based on Program PLOT. FOR. Dec 1987 by M.A. Morgan 
C 
C INPUT DATA FORMAT 
C 
C TITLE - 64 Space Title 
C NPTS - # Data Points 
C XMIN - Real Min XK value 
C XMAX - Real Max K value 
C F(N) - Input Data Array 
C IHC - Hardcopy (1=Yes) 
C 


CHARACTER*1 YN,DUM 
CHARACTER*64 TITLE,XS,ID,YS,RELAT 
REAL X(1025),F( 1025) ,XMIN,XMAX,FFF( 1025) ,XMINN, XMAXX, L9 
INTEGER*2 N,JROW,JCOL, ISYMBL,ITYPE,CYAN, WHITE, YELLOW 
IHC=0 
WRITE(*,*) ' WOULD YOU LIKE TO HAVE A HARD COPY: ' 
READ(*,100) YN 
IF(YN. EQ. 'Y'. OR. YN. EQ. 'y') IHC=1 
WHITE=7 
CYAN=11 
DO 1111 I=1,NPTS 
F(1)=FFF(T) 
1111 CONTINUE © 
XMIN=XMINN 
XMAX=XMAXX 
YELLOW=14 
=NPTS 
DX=( XMAX-XMIN)/(NPTS-1. 0) 
FMIN=0. 0 
FMAX=0. 0 
DO 22 K=1,NPTS 
X(K)=XMIN+(K-1. 0)*DX 
IF(F(K). LT. FMIN) FMIN=F(K) 
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IFCFC(K).GT. FMAX) FMAX=F(K) 
ZZ CONTINUE 
IFC FMIN. GT. 0.0) FMIN=0. 0 
IF( FMAX. LT. 0.0) FMAX=0. 0 
Computing Scale Factors for Vertical Axis 
ABSMIN=ABS(FMIN) 
ABSMAX=ABS ( FMAX) 
YMAX=AMAX1( ABSMIN, ABSMAX) 
NSCL=INT( LOG10( YMAX) ) 

IF (YMAX. LT. 1.0) NSCL=NSCL-1 
YSCL=10. **NSCL 
FMIN=FMIN/YSCL 
FMAX=FMAX/YSCL 

ABSMIN=ABSMIN/YSCL 
ABSMAX=ABSMAX/YSCL 
DO 33 K=1,NPTS 


33 F(K)=F(K)/YSCL 
YMIN=0. 0 
IF(FMIN. EQ.0.0) GO TO 37 
35 YMIN=YMIN+0. 5 


IFCABSMIN. GT. YMIN) GO TO 35 
YMIN=YMIN*FMIN/ABSMIN 


oe CONTINUE 
YMAX=0. 0 
IFC FMAX. EQ.0.0) GO TO 41 
39 YMAX=YMAX+0. 5 


IF( ABSMAX. GT. YMAX) GO TO 39 
YMAX=YMAX**FMAX / ABSMAX 
41 CONTINUE 
CALL QBEEP 
CALL QSMODE(2) 
WRITE(*,*) ‘Press RET to View Plot...RET Again to Clear Screen’ 
IFC IHC. NE. 1) GO TO 42 
WRITE(*,*) 'A Hardcopy Will Be Made...Check that Printer is ON’ 
42 CONTINUE 
READ(*,100) DUM 
Calling GRAFMATIC Routines and Plotting Solid Line Graph 
ITYPE=1 
ISYMBL=-2 
CALL QSMODE( 16) 
CALL QPTXT(64,TITLE, YELLOW, 16,23) 
CALL QPTXT(64,ID,CYAN,16,19) 
CALL QPTXT(64,XS ,CYAN, 16,20) 
CALL QPTXT(64,YS ,CYAN, 16,21) 
RELAT=' RELATIVE SMOOTH =' 
CALL QPTXT( 64 ,RELAT,CYAN,16,18) 
CALL QQNPUT( 265 ,253,L9,3) 
CALL QPLOT( 110 ,540,20,240, XMIN, XMAX , YMIN, YMAX,0. ,0. , 0 eee 
CALL QSETUP(0,CYAN, ISYMBL, CYAN) 
XMAJOR=XMAX/5. 0 
CALL QXAXISCXMIN, XMAX ,XMAJOR,1,1,2) 
YMAJOR=1. 0 
CALL QYAXIS(C YMIN, YMAX, YMAJOR,1,1,2) 
JROW=14+( ABSCYMIN) /CABS( YMAX )+ABS( YMIN) ) )*220 
JCOL=80-430*XMIN/( XMAX -XMIN) 
CALL QGTXT(3,'0.0' , WHITE, JCOL, JROW,0) 
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44 


100 
150 


GALINOPIX101, s_, YELLOW ,5,16) 
Sl Ormat(1, c ,YELLOW,5,15) 
CALL QPTXT(1,'a' ,YELLOW,5,14) 
GAbL OPGxr(1, 1° , YELLOWS, 13) 
GAO RExrC1, e| ,YELLOW.5 12) 
ORE QPTXT( 1, 'X' , YELLOW,5, 10) 
CALL QPTXT(2,'10° , YELLOW,4,8) 
CALL QPTXT(1,'*' , YELLOW,5,7) 
CALL OPmxt(1,° *”,YEDLOW,5 ,6) 
CALL QCMOV(0,6) 
WRITE(*,150) NSCL 
CALL QTABLCITYPE,N,X,F) 
READ(*,100) DUM 
IF(IHC.NE.1) GO TO 44 
CALL PRTSC 
CONTINUE 
CALL QSMODE( 2) 
FORMAT(A) 
FORMAT( 4X,I3) 
RETURN 
END 


SUBROUTINE PLO77(Y1,Y5,T2,NUMPT, ID1) 


Cra tedetedsdevede tess eve se Tete Lede Tee Ke Fe Fe Te Te TeVET TENT TENE HATER ERE REE ER 


C 
C 


Y2 = Y5 = BACKGROUND DATA ARRAY * 
yal = SCATTERED DATA ARRAY * 


Crete dededvedede de evede FeWe ye Te Fe Fe Fe Fe Te Tee Ts PETC TEN Te TCT TENE NE NIE EERE ERT ER EEE 


Ze 


INTEGER NUMPT , DATAB ,SUBAV,F1,NSHIF,SET 

REAL Y1( 1024) ,Y2( 1024) ,¥3( 1024) ,Y4( 1024) , Y5( 1024) 

REAL XMIN,XMAX 

CHARACTER*16 NAME1,NAME2 ,NAME3 ,QN 

CHARACTER*1 YN 

CHARACTER*64 YS,XS,ID,TITLE,DATE 

CHARACTER*64 ID1,TYPE1 

CHARACTER*64 ID2,TYPE2 

CHARACTER*64 ID3,TYPE3 

REAL TIME,T2 

WRITE(* ,*) " dovledevedevevevedevedededete devedete de Leeeke Te E TRL NM 
WRITE(*,*) ‘SHIFTING SCATTERED WAVEFORM & SUBTRACTING BACKGROUND' 
WRITEC* ,**) ' Seveverdedededevedededededeieledededeseledesetesedesese ete sete deed ede Fee Ge PERRET ESI 
XMIN=0. 0 

XMAX=T2 

DO 22 K=1,NUMPT 

Y2(K)=Y5(K) 

SET=0 


CreredtedededeteKR RRR TCR RTE RL LR RETR REEVE RERERR GE 


33 CALL SHIFTCTITLE,YS,XS,ID) 
Crete tededecede dese cece dese se Feed ede TOT ED CTT. TFC ETE TERRELL TR LR ETRE 


CALL PLO8( TITLE, NUMPT, XMIN, XMAX,Y1,XS,YS,Y2, ID) 


Credededededesede sede Teche de he Tete eRe TRL ERATE LTTE EREERELTETERR EG 


C 


SHIFT THE CURVE 


Cretetere dese dete se fered Fee RT RETR PRR REECE RRER REE RG 


CALL NSHIFT(NSHIF) 

IF(NSHIF. EQ.0) GO TO 40 

CALL GOOD(Y1,Y2,NSHIF ,NUMPT, SET) 
GO TO 33 
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40 CONTINUE 
CreVeHe Fea Severs Kete Me Fe R TCT TK IK IKI KAERERERERKE 


€ SUBTRACT BACKGROUND FROM THE WAVEFORM ix 

CAH Te Heese Fee Tee Tout Fese Fe Tee Te Ne Fe He TE HN RT Tee TK RTI RAIREEEREEREER 
DO 44 I=1,NUMPT 

44 Y3(1I)=¥1(1)-Y2(1) 
CALL SUBTR(TITLE,XS,YS) 
CALL P11(TITLE,NUMPT , XMIN,XMAX,Y3,XS,1D1,YS) 
WRITE(*,*) ‘DO YOU WANT TO RE-SHIFT AND SUBTRACT ? (Y/N): 
READ(**,300) YN 
IF(YN. EQ. 'y'. OR. YN. EQ. 'Y') GO TO 33 
DO 77 IJ=1,NUMPT 
77) =—-Y1C1)=Y3¢1) 
300 FORMAT(A) 

RETURN 
END 


C aeveteteseve dee te He LTE TET LOTS SE TETS FTE TEESE Le TE SE SO TE LED VE TE EIS TTC SE TCT SEE HE 


SUBROUTINE NSHIFTCNSHIF) 

INTEGER NSHIF 

WRITE(*,*) UV Sedede dered de dese dese dese sede IT NIRA | 
WRITEC* ,**) 

WRITE(*,*) ‘SCATTERED WAVEFORM SHIFTING:  ' 

WRITE(* ,* ) 

WRITE(*,*) ‘+n shifts scat waveform to the RIGHT n-points' 
WRITE(*,*) '-n shifts scat waveform to the LEFT n-points' 
WRITE(*,*) ' O DISPLAYS subtracted waveform 

WRITE(*,*) 

WRITE(*,*) 'ENTER SHIFT (+/-INTEGER):  ' 

READ(*,*) NSHIF 

RETURN 

END 


Ceevede deseo cede secede Te Fede Tee Te Te BETS POTEET DELETE TEETER TERETE TETOLEE RTGS 


SUBROUTINE SUBTR(TITLE,XS,YS) 
CHARACTER*64 TITLE,XS,1D,YS 
TITLE='SUBTRACTED WAVEFORM' 
XS='X-AXIS TIME IN NS' 
YS='Y-AXIS MAGNITUDE IN VOLT' 
RETURN 

END 


Cvededededededede ve de tess Fede Fes e Te Fe TEE TELE TE TER TELE LICE TENE LL TEKST REEL 


SUBROUTINE SHIFT( TITLE, YS,XS,ID) 
CHARACTER*64 TITLE,YS,XS,ID 
TITLE='TIME SERIES PLOT' 
YS='BLUE-TARGET' 
ID='GREEN-BACKGROUND' 
XS='SHIFT-TARGET ONLY’ 

RETURN 

END 


C vedo devededede tod dete de cede Tove ese Fe ToS PTE OEE CTE TELE TE LES TELE. SENT TER E S 


SUBROUTINE GOODCY1,Y2,NSHIF,NUMPT,SET) 


GC ve 
C PROGRAM TO SHIFT Y2 CURVE * 
C 

C = SHE he YL. TO) THE LEFT * 
C + SHIFT Y1 TO THE RIGHT 7 
C 

C er IS THE SHIFTING CURVE ¥ 
C 1G: IS THE BACKGROUND CURVE * 
C NSHIF IS THE NUMBER OF POINT TO SHIFT * 
C SET IS THE SUM OF THE POINT TO SHIFT * 
C IN ORDER TO ZERO PADDING * 
C NUMPT IS THE NUMBER OF POINT ¥ 
C v 
Crteasedeaden delisted tics texte Fetes Ha Hee He Fe Fe Fe HF Te He eH HH eK KK Te HH 


REAL Y1( 1024) ,Y2( 1024) ,Y3(1024) 
INTEGER NSHIF,NUMPT,SET 
5 IF( NSHIF.LT.O ) THEN 
NSHIF=- 1*(NSHIF) 
DO 20 I=NSHIFT,NUMPT-NSHIF-SET 
Y3(1)=Y1(I+NSHIF) 
20 CONTINUE 
DO 10 I=NUMPT-NSHIF-SET,NUMPT 
¥2(1)=0. 0 
Y3(1)=0. 0 
10 CONTINUE 
DO 15 I=1,NSHIF+SET+1 
VOCUS. 0 
WACO) 0 
15 CONTINUE 
ELSE IF( NSHIF.GT.O ) THEN 
N=NSHIF+SET+1 
DO 40 I=N,NUMPT 
Y3(1)=Y1(1-NSHIF) 
40 CONTINUE 


30 CONTINUE 
DO 35 I=NUMPT-N 
Y2(1)=0.0 
Y3(1)=0.0 
35 CONTINUE 
ELSE 
END IF 
SET=SET+NSHIF 
DO 50 I=1,NUMPT 
Y1(1)=Y3(1) 
50 CONTINUE 


SHIF-SET ,NUMPT 


RETURN 

END 

SUBROUTINE OUTCNAME ,ID,TYPE,DATE,F1,SUBAV,T2,DATAB,NUMPT,Y) 
C 
C ID $ = CHARACTER $ TARGET IDENTIFICATION 


65 


SP ep ep ep | ap (eye ap aap, 


10 
20 
Zu 
Qe 
30 


TYPE $ = TYPE OF WAVEFORM 

DATE $ = DATE OF EXPERIMENT TAKE PLACE 
Fl = DUMMY STRING 

SUBAV = SUB AVERAGE 

NUMPT | = NUMPT OF POINT 

DATAB = DATA BLOCK 

Y = Y ARRAY 

T2 = TIME WINDOW 


CHARACTER*64 ID,TYPE,DATE,F1 
CHARACTER*16 NAME 
REAL Y( 1024) ,T2 
INTEGER SUBAV,NUMPT,DATAB 
OPEN( UNIT=10 , FILE=NAME , STATUS=' UNKNOWN' ) 
WRITE(10,*) ID 
WRITE(10,*) TYPE 
WRITE(10,*) DATE 
WRITE(10,*) Fl 
WRITE( 10,22) SUBAV 
WRITE( 10,21) T2 
WRITE( 10,22) DATAB 
WRITE( 10,22) NUMPT 
DO 10 I=1,NUMPT 
WRITE(10,20) Y(I) 
CONTINUE 
FORMAT(E12. 6) 
FORMAT(F5. 2) 
FORMAT(I5) 
FORMAT( A) 
CLOSE(10) 
RETURN 
END 
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APPENDIX B. 


THE ACQUISITION QUICK BASIC PROGRAM 


The following computer program acquires the measurement data from the HP 
54120T as described in chapter 5. The program 1s written in QUICK BASIC 


1000 
1020 
1040 
1060 
1080 
1240 
Te21 
22 
1324 
TZ 5 
20 
roy 
IL Sy/Ate: 
1329 
1330 
1341 
1342 
1343 
1344 
1345 
1346 
1347 
1360 
1380 
1400 
1420 
1440 
1460 
1480 
1500 
1520 
1540 
1560 
1580 
1600 
1620 
1640 
1660 
1680 
1700 
1720 
1740 
1760 
1780 
1800 
1820 
1840 


Fetedededesesede de dede Feces Reve RRR RIITRRETTRT RERREERERER EE 


MODIFY 4/22/88 
WRITTEN BY LT. JG SOONPUEN SOMAPEE 
ACQUISITION 


REM 
REM * 
REM * 
REM * 
REM * 
REM vc Je 
REM Sevevedededtedesvevedevededeklecesedescdesedescdesedesesdcksesdedescddekkkvctkedeve 
DS=DATES 
DUMMYS=TIMES: PRINT ''DATE..... 
COLOR 7,1 
PRINT "WELCOME TO ACQUISITION PROGRAM" 
PRINT ''THERE ARE 3 TYPE OF MEASUREMENT" 
PRINT “1. BACKGROUND...... WAVETYPE'' 
PRINT "2. CALIBRATE....... WAVETY PE" 
PRINT ‘3. TARGET.........- WAVETYPE"' 
PRINT '' PRESS THE FOLLOWING KEY " 
PRINT " T ( USED DATA "LAST RESPONSE).............. 4 
PRINT " D (lS) DEBAULM mms ee... - 
VS=INKEYS: IF VS=""' THEN 1344 
IF VS="L" OR VS ="'1'' THEN 1960 
IF vS="d" OR VS ="D" THEN 1360 
GOTO 1344 
REM Sevevevesedesesotedestesedisesededsededesedetevedetesestedeietetekdlesketetedseeskses: 
REM * DEFAULT ve 
REM Fi RRR OIRO iis Ria fant telnet tale toie port airfise (sr f ur ive ish t orice ise rie heehee berl or hse fori 
OPEN "DATA. DAT'' FOR OUTPUT AS #1 
P1=512 
PRINT #1,P1 
N=1024 
PRINT#1,N 
ND=1 
PRINT##1 ,ND 
TS="4 inch Thin wire" 
PRINT#1,TS 
DS=DATES 
PRINT#1,DS 
T2=1. 563E-10 
PRINT#1 ,T2 
WS="Background Wave Type" 
PRINT##1,WS 
F1=20 
PRINT#1,F1 
Sia 
PRINT#1,S1 
G1=Sl 
DSES="Backg' 
PRINT##1,DSES 


+  & of 


";DS: PRINT ''TIME..... "'; DUMMYS 


67 


1860 
1880 
1900 
1920 
1940 
1960 
1980 
2000 
2020 
2040 
2060 
2080 
2100 
2120 
2140 
Z160) 
2180 
7200 
2220 
2240 
2260 
2280 
2300 
Za20 
2340 
2360 
2380 
2400 
2420 
2440 
2460 
2480 


DUMMY S=TIME 


: 


PRINT#1,DUMMYS 


DELAY 1=38 


PRINT#1, DELAY1 


CLOSE#1 


REM Sorcverkvokokeaksakovekseskseakskoe koko skveseokedeokokakscdeseveskokakeseokseseseskseste se 


REM 
REM 


Optionally include the following declarations 
in your program. 


REM They provide appropriate mnemonics by which 


REM 
REM 


to reference commonly used values. 
Some mnemonics (GET%, ERR%, 


REM END%, ATN%) are preceded by "B" in order to 
REM distinguish them from 


REM 
REM 
REM 
UNLZ 
UNT% 
GTL% 
SDC% 
PPC% 
BGET% 
eA 


foo UW 


&H3F 
&HSF 
&H1 
&H4 
&H5 
= &H8 
&H9 


LLO% 
DCL% 
PPU% 
SPE% 
SPD% 
Per 
PPD% 
REM 


nti id du neu 


&H11 
&H14 
&H15 
&H18 
&H19 
&H60 
&H70 


GPIB Commands 
' GPIB 


t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 
t 


BASICA keywords. 


unlisten command 
untalk command 

go to local 

selected device clear 
parallel poll configure 
group execute trigger 
take control 


GPIB 
GPIB 
GPIB 
GPIB 
Gras 
GPIB 


GPIB 
GPIB 
GPIB 


GPIB- 


GPIB 
GPIB 
GPIB 


local lock out 
device clear 
ppoll unconfigure 


serial poll enable 


serial poll disable 
parallel poll enable 
parallel poll disable 


REM GPIB status bit vector 

REM global variable IBSTA% and wait mask 
Error detected 
Timeout 

EOI or EOS detected 
SRQ detected by CIC 
Device needs service 
I/O completed 

Local lockout state 
Remote state 
Controller-In-Charge 
Attention asserted 
Talker active 
Listener active 
Device trigger state 
Device clear state 


t 
t 
t 
{ 
t 
t 
{ 
t 
t 
t 
t 
t 
t 


Reh Error messages returned in global variable IBERR% 

DOS error 

Function requires GPIB-PC to be CIC 
Write function detected no Listeners 
Interface board not addressed correctly 
Invalid argument to function call 
Function requires GPIB-PC to be SAC 

I/O operation aborted 


BERR% = &H8000 
TIMO% = &H4000 
BEND% = &H2000 
SRQI% = &H1000 
RQS% = &Hco00 
CMPL% = &H100 
LOK% = &H80 
REM% = &H40 
CICK = &H20 
BAIN‘ = &H10 
TACS% = &H8 
LACS% = &H4 
DTAS% = &H2 
DCAS% = &H1 
REM 
EDVR&Z = 0 
ECICZ = 1 
ENOLZ = 2 
EADR%Z = 3 
EARG%Z = 4 
ESAC% = 5 
EABO% = 6 


t 
t 
t 
t 
t 
t 
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2980 
3000 
3020 
3040 
3060 
3080 
3100 
I 120 
3140 
3160 
3180 
3200 
3220 
3240 
3260 
3280 
3300 
3320 
3340 
3360 
3380 
3400 
3420 
3440 
3460 
3480 
3500 
3520 
3540 
3560 
3580 
3600 
3620 
3640 
3660 
3680 
3700 
3720 
3740 
3760 
3780 
3800 
3820 
3840 
3860 
3880 
3900 
3920 
3940 
3960 
3980 
4000 
4020 
4040 
4060 
4080 


ENEB% 
EOIP% 
REM 

ECAP% 
EFSO% 
EBUS% 
ESTB% 
ESRQ% 
REM 

REM EOS 


7 t 
10 


11 
IZ 
14 
15 
16 


mode bits 


BIN% = &H1000 


XEOS% = 
REOS% = 
REM 


&H800 
&H400 


Non-existent interface board 

I/O operation started before previous 
operation completed 

No capability for operation 

File system operation error 

Command error during device call 
Serial poll status byte lost 

SRQ remains asserted 


Eight bit compare 
Send EOI with EOS byte 
Terminate read on EOS 


REM Timeout values and meanings 


TNONE% = O Infinite timeout (disabled) 
T10US% = 1 ' Timeout of 10 us (ideal) 
T30US% = 2 ' Timeout of 30 us (ideal) 
T100US% = 3 ' Timeout of 100 us (ideal) 
Re0ous, = 4 ' Timeout of 300 us (ideal) 
T1MS% = 5 ' Timeout of 1 ms (ideal) 
T3MS% = 6 ' Timeout of 3 ms (ideal) 
T10MS% = 7 ' Timeout of 10 ms (ideal) 
T30MS% = 8 ' Timeout of 30 ms (ideal) 
T100MS% = 9 ' Timeout of 100 ms (ideal) 
T300MS% = 10 ' Timeout of 300 ms (ideal) 
hs. = dal PTameout of 1 s (ideal) 
T385% = 12 ' Timeout of 3 s (ideal) 
T108S% = 13 ' Timeout of 10 s (ideal) 
T305% =" 14 ' Timeout of 30 s (ideal) 
T100S% = 15 ' Timeout of 100 s (ideal) 
T3008S% = 16 ' Timeout of 300 s (ideal) 
TIGOOS,% = 17 ' Timeout of 1000 s (maximum) 
REM 

REM Miscellaneous 

S% = &H8 Parallel Poll sense bit 
LF% = &HA ' Line feed character 

REM 


REM Application program variables passed to 
REM GPIB functions 


REM 

CMDS = SPACES(10) ' command buffer 

RDS = SPACES( 255) ' read data buffer 

WRTS = SPACES(255) ' write data buffer 

BNAMES$ = SPACES(7) ' board name buffer 

BDNAMES = SPACES(7) ' board or device name buffer 
FLNAMES = SPACES(50) ' file name buffer 


REM Setedevededededededenodedesedededededetededededetededekdedtesedetetededetetetede 


CLS 
DIM W( 1500) ,RRR( 1500) 
DIM RDSS$( 1024) 


SCREEN 0 

COMOR® 7,1,3 

OPEN "C: DATA. DATA" FOR INPUT AS #1 

BEEe: CLS 

REM  tiesdedededevedededededeveveveteveveveds Fetededese dedede detededeve de de ted sede he 
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REM * BEGIN MAIN PROGRAM * 

REM wacdededededkacdededevoue ded devedenke deve deve dere de de deve dese deve Hevea Hex sex 

LOCATE 5,8: PRINT’ DATA ACQUISITION 

REM Weede dese devede de devede eve de Fe vee He ve Teds ae ae ak ae He Ve Fe He Te Fe Fe He Ve Te se Tose He He Te se ae Fe Te Fe Te Fe ee We He He He te Heese Kee 
REM * TS = TARGET IDENTIFICATION * 
REM * WS = WAVEFORM TYPE * 
REM * DS = DATE * 
REM * N = NUMBER OF SUB-AVERAGE / WAVEFORM OR COUNT * 
REM * T2 = DPO SAMPLING INTERVAL (nsec) * 
REM * S1 = DPO MAXIMUM DPO VERTICAL SCALE(M) * 
REM * ND = NUMBER OF DATA BLOCKS * 
REM * Pl = NUMBER OF POINT / WAVEFORM * 
REM * Fl = DPO TIME WINDOW * 
REM * DSES = FILE NAME OF THE DATA OUTPUT FILE 7 
REM * DEVS = NAME OF THE HPDEV IN THE CONFIGURATION FILE* 
REM * DSO% = THE NUMBER ASSOCIATED WITH THE DEVICE 

REM * VEWS = STRINGS VARIABLE FOR SENDING TO HP54120 * 
REM * DUMMYS = TIME SET * 
REM  xaedededcdededededededede deve teas de deve deve ve ve de de ve de Fede He ve Weve Weve Tee Ne Tove de Teves He Fee VEN Ie Fe Tere We HA Ve 
REM 

INPUT#1,P1 

PRINT "NUMBER OF POUNG.).. 0 eer ". Pal 

INPUT#1,N 

PRINT "NUMBER OF SUB-AVERAGE(count)....... ".N 

INPUT#1 ,ND 

PRINT "NUMBER OF DATA BLOCK............... ". ND 

INPUT#1,TS$ 

PRINT “TARGET ID.... ee a: @ sce) eee ":TS 

INPUT#1,D$ 

PRINT “DATE....:22e:. . (eee "3 DS 

INE Oils a2 

PRINT "DPO SAMPLING INTERVAL.............. et We 

INPUT#1, WS 

PRINT? ss WAVEFORMes GYPE.. >... center see "s Ws 

INPUT#1,F1 

PRINT "“DPO TIMEWINDOW( NSEC )......... ". BY 

INPUT#1, $1 

PRINT “DPO VERTICAL SCALE@MO mee... ee ncn 

Gi=sl 

INPUL ee Dors 

PRINT "DATA FILE NAME OUTPUT.............. "; DSES 

INPUT#1 , DUMMYS 

PRINT “AUTO TIME QM. 2... fee... .. eee ''; DUMMYS 
INPUT#1 , DELAY1 

PRINTS DELAY <C NSEC ) sige cee ". DELAY1 

CLOSE#1 

PRINT oo ccc oo 0 me ace ec Ue ce nce 
COLOR 7,1: PRINT '' PRESS THE FOLLOWING KEY " 

PRINT " Y ( CHANGE THE DATA. DAT ).-...«..- Seg. - Jee a 
PRINT '' N ( NO CHANGE )...... 2.0. % 0000 Smee en oe " 


VS=INKEYS: IF VS=""' THEN 5140 


IF VS="Y" OR VS 
IF VS="'N" OR VS 


"y" THEN 11260 
"n" THEN 5220 


GOTO 5100 
DEVS="HPDEV' 


REM 


seslevevevesededevedtedevededevededesesetovetedesecesedesvedkidsedesedledesetedkdevek 
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w 


5260 
5280 
5300 
Do 20 
5340 
5360 
5380 
5400 
5420 
5440 
5460 
5480 
5500 
5520 
5540 
5560 
5580 
5600 
5620 
5640 
5660 
5680 
5700 
5720 
5740 
5760 
5780 
5800 
5820 
5840 
5860 
5880 
5900 
20 
5940 
519,60 
5980 
6000 
6020 
6040 
6060 
6080 
6100 
6120 
6140 
6160 
6180 
6200 
6220 
6240 
6260 
6280 
6300 
6320 
6340 
6360 


** OPEN DEVICE AND RETURN THE UNIT ve 
* DESCRIPTOR ASSOCIATED WITH THE GIVEN * 

REM * NAME * 

REM sevevescvescsededededevededeseseactdededevevedetededesvetededede vee de desedede vse ve 

CALL IBFIND(C DEVS , DSO%) 

IF DSO% < 0 THEN PRINT "ERROR IN IBFIND" 

REM  weverededededededededededevedede se vede te ve cede Tove vee Fe Fe Fe Fe Fee Hee Hee eI He 


REM * CLEAR SPECIFIED DEVICE * 
REM atdededededodedesdedesededeaededeuededene He ve reve Fever ve Fede ve Fete ve He seve He ve Fe Fe 


CALL IBCLR(DSO%) 


REM 
REM 


REM xaevevededededededesedede dedede dese deve dee Fede seve Fe Te Fe He Hee se Fee Fe He Fe Fe HE 
GOSUB 6080 

CLS: BEEP 

PRINT WS+'' DATA STORAGE COMPLETE" 


PRINT '' MEASURE OTHER WAVEFORM TYPE(Y/N) " 
QS=INKEYS: IF QS=" " THEN GOTO 5560 
IF QS="N" OR QS="n" THEN GOTO 5640 


IF QS="Y" OR QS="y"' THEN GOTO 4040 

GOTO 5560 

SOP 

RE Mivrereavsedevedevevede deve deve dove se vevede ve ve deeded te Fe ve Fede Tose He Fee Te Fe He He 
REM* SUBPROGRAM * 
RE M* * 
REM Sedeakvededededesededdodoeivededededescdedesdesedesedesesedededescicdedededesedede 
REM * TEXT REFERENCE ¥ 
REM * "GURU II GPIB USER'S RESOURCE" * 
REM * "UTILITY FOR IBM PC" * 
REM * "TEKTRONIX 1986" * 
REM _ waeatdedededededsdedcdede dedeat dese de sede de vede Fee se te ac de severe ve FoI Fede ve seve 
REM x veatvedeaess dese devedoceledevevedesk deaede seus Has ve deat sede x ve Fe Fee He ve Te ve ve 
REM * WRITE DATA FROM STRING SUBROUTINE ¥ 
REM _ tevededsasdedevedovevevevevesvaedeveve seve cess seve de dese seve ve ve eve Fe ve ve Fe Fe Ie 7 
CALL IBWRT(DS0%, VEWS) 

RETURN 

REM  wWaerededededededesestcaede doves ae ve devededevedede se deve Fede FeIe Fe Fede Te Fe Fe Fe He Fe 
REM * TEXT REFERENCE 7% 
REM * "GETTING STARTED GUIDE" ¥ 
REM * "HP 54120T DIGITIZING " ¥ 
REM * "OSCILLOSCOPE" * 
REM * "HP COMPANY 1987" * 
REM  xevededededcvedededescvesedeveve sede se vewevewe ve seve ve te se seve se ae Fe ve vege IF Ie Fe Ve 
REM wevededevevededevedevevedese reds Tete Te Bede TEER. PETE TERETE DE TEE PETC SCTE TET 
REM * SETUP THE DISPLAY * 
REM wetevedededevededededevese deve cede ceded eae Fe TEEN TCT RR LTCR REET 
VEWS="'*RST'! 


GOSUB 5840 

REM xredededededesededededededededededededeaede de dese seven dene deed sede de desede Ke 
REM * CRT FORMAT 1 PROVIDE ONE DISPLAY AREA * 
REM * AND USE 8 DIVISIONS FOR THE FULL ¥ 
REM * SCALE RANGE ¥c 
REM x xeaesevestrededededededededede dese rededede dese re deve deve de sevede dete Heve Ke vedere 
VEWS="': DISPLAY: GRATICULE FRAME; FORMAT 1” 


GOSUB 5840 
REM 3 wrvevededededededededevedededeseKe te RRR RR ERRREREECRER I 


REM * TURN TDR STEP GENERATOR ON ¥ 
REM Jededevevedvevedevevevevedevevevedevevevededesvevedetesesesedededededesekkictesede 


7] 


6380 
6400 
6420 
6440 
6460 
6480 
6500 
6520 
6540 
6560 
6580 
6600 
6620 
6640 
6660 
6680 
6700 
6720 
6740 
6760 
6780 
6840 
6860 
6880 
7020 
7040 
7060 
7080 
7100 
720 
7140 
7160 
7180 
7200 
7220 
7240 
7260 
7280 
7300 
7320 
7340 
7360 
7280 
7400 
7420 
7440 
7460 
7480 
7500 
7520 
7540 
7560 
7580 
7600 
7620 
7640 


VEWS="': NETWORK: REFLECTION: PRESET" 

GOSUB 5840 

REM ttededkevedesededeoesede eae dese ds deve de ak eve de Towed seve dees Fe He He eH He Me Fe HK 
REM * PAGE 6-13 TURN ON CHANNEL 4& * 
REM  vevedevesedecesededededs vedeuevededeakakakede sted ak ak seve ae aeak Hedk seeded ve ses oe 
VEWS=": VIEW CHANNEL4" 

GOSUB 5840 

REM  tevedededededededsvbedeakvkevodesonese se douse sede deve dee dee ve Me se Tes We ve He He Tee 
REM * PAGE 6-4 TURN OFF CHANNELI * 
REM xtedededsdededsse dese dededevedc deve deve ve de dese te Heke dete Teak ek HHH ke He ee HK 


VEWS="': BLANK CHANNEL1" 

GOSUB 5840 

REM  teredededecededededededesededk ved aedevedene Heed Wea Vek HTK Kee Kee eH 
REM * SET UP CHANNEL RANGE AND OFFSET % 
REM * PAGE 17-6 * 
REM * TO BE 20 NANOSECONDS * 
REM tered dedevededesede devon deaedevesede souk aeakat Weak dere Here Vere HH Kee Nk He 
VEWS="': TIMEBASE: RANGE 20 NS" 


GOSUB 5840 


REM PRINT WM eeseatsestsvokakslecskvevedeskseskvedeaesedoskedeskieseskakokacdedesesededededcdedeskedesede aesedeaeskes |! 
REM PRINT "SET THE VERTICAL SCALE RANGE FROM 8 TO 640 MV 


VERT=S1 
VEWS="': CHANNEL4: RANGE ''+STRS(VERT)+'M' 
GOSUB 5840 


REM wrededcdedercdkesevevedede dese de vede sede Weve de eee NTE ee Ie VEIN He HEH He 
REM * TURN ON THE GRATITUDE * 
REM  tededededodevevestdededededeaedede de dese dee ae ods seve Fete de Ne He Vee dete HIE He 
VEWS="': DISPLAY: GRATICULE GRID" 

GOSUB 5840 

VEWS="': ACQ: BAND LOW" 

GOSUB 5840 

REM xtvedededevedesedeskdovedededevedede deve ve dene dee de Fe ve Tove He He eI Hee He HH HK 
REM * — AVERAGING ON ¥ 
REM : ste ve euteate ese a'e s'e te we wes s'e ede destes) s'e reve s‘c steves! ale ewe ste Sele seste seslesleclesesevovc we 
ave" -acourre. ae AVERAGE" 

GOSUB 5840 

REM  wredsedodedtobv deck vedededksede ded dese ve dove se ae deve ve dese Weve He Tete Weve Heese H 
Jee) oil as ASKING FOR BANDWIDTH ve 
REM  xveakeaesestcedsscakiede dese de devas ve dk ve de vie ak deals He de ake ae ake ae He We Hee Heke ve 
PRINT tee a Sleleclewtontoateatectoctonterte cteatoatonle leatealoctentoctentoctecleatontc ate desloge slesledesie se t? 
PRINT "THE DEFAULT BANDWIDTH IS 12. 4GHZ ? 
PRINT "WOULD YOU LIKE TO CHANGE TO 20 GHZ (Y/N)" 


QS=INKEYS: IF Q$=" " THEN GOTO 7380 
IF QS="N" OR QS="'n'' THEN GOTO 7500 
IF Qs="Y" OR QS$="y'"' THEN GOTO 7460 
GOTO 7380 

VEWS=": ACQ: BAND HIGH" 

GOSUB 5840 

RDS=SPACES( 255) 


VEWS="ACQUIRE: COUNT ''+STRS(N) 

GOSUB 5840 

VEWS="'; ACQUIRE: POINTS ''+STRS(P1) 

GOSUB 5840 

REM xtededededsdecedededededede detec de ded Ke RRR REECE 1X 
REM * SbL TIME SDEEAY * 
REM * Py 


qe 


7660 
7680 
7700 
7720 
7740 
7760 
7780 
7800 
7820 
7840 
7860 
7880 
7900 
7920 
7940 
7960 
7980 
8000 
8020 
8040 
8320 
8340 
8360 
8380 
8400 
8420 
8560 
8600 
8620 
8640 
8660 
8680 
8700 
8720 
8740 
8760 
8780 
8800 
8820 
8840 
8860 
8880 
8900 
8920 
8940 
8960 
8980 
9000 
9020 
9040 
9060 
9080 
9100 
9120 
9140 
9160 


REM Mires tes Kes sewese te sess sesese te sees dese seve se see ye se se Fe Fe Ne eM 


REM 

REM 

REM 
DELAY=DELAY1 
DUMS$=STRS ( DELAY) 
VEWS=": TIM: DEL ‘'+DUMS+"NS" 

GOSUB 5840 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM wervdededededededededededede dese sede de desede de deve Heve He Fe Fe Hee HK Ie Ie He HH 

REM HP 54120T SENT THE LAST BYTE RESPONSE "LF'' ASII DECIMAL 10 
VEWS="': SYSTEM: HEADER OFF;:EOI ON" 


GOSUB 5840 
PRINT "Vaevedeovcvdescoereaeaeaeobdedese dese dedeak sede de sede seve de ak dese de de sede Weve desea ese se Wes Ie ese 


REM eres se se Vee Ve Ve TE VES eNeIETE VEIE TE EI IE TO NE IO TE IE IT TENE TI IE IT IT IIIT 


PRINT "I AM GOING TO DIGITIZE CHANNEL( Y/N)" 
QS=INKEYS: IF QS=" " THEN GOTO 8360 

IF QS="Y" OR QS="y" THEN GOTO 8420 

GOTO 8360 

PRINT TIMES 


REM wreverevevevevedestedeuedesesc sede seak deseve deve desk ac ae dese sie ae ve oe ve ae ve desk deat 


REM * GET THE NUMBER AND STORE IN THE * 
REM * FILE * 
REM wrvtedededssededese deve tesedevedede dese dese seve seve Wee ese ve se se se ve He se He Fes ve 
DSEWS="'C: C" 

FOR I=1 TO ND 


IF I<10 THEN DSFILES=DSEWS+RIGHTS(STRS(1I),1)+". WFM" 
IF I>9 THEN DSFILES=DSEWS+RIGHTS(STRS(1),2)+'". WFM' 
REM 

REM wtevededevedevevededevedede dete sede Kec Te Pe T CLEP TE TE TCR TEESE TEPER TEVETEIE 

REM * DIGITIZE CHANNEL 4 * 

REM wvevededevodevovedevovededevede decade dese tect Feces Be ICTE RTT REIKI 


VEWS="': DIGITIZE CHANNEL4" 


REM 

GOSUB 5840 

VEWS="': WAVEFORM: SOURCE WMEMORY4; FORMAT ASCII" 
GOSUB 5840 

VEWS=": WAVEFORM: DATA?" 

GOSUB 5840 

REM 

CALL IBRDF(DS0%,DSFILES ) 

CLOSE#1 

NEXT I 

REM TePeTH RK RRR TERRI LET VER ER LIRR IRR RRR 
REM * FINISH DIGITIZING * 
REM  wrevedevedevedededededededese ds devede ce Fede RRR RR RR RICK 
PRINT TIMES 


DSFILES=""HEADER" 
VEWS=": WAV: PRE?" 
GOSUB 5840 
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CALL IBRDF(DSO%, DSFILES ) 

PRINT "HEADER COMPLETE” 

GOSUB 10200 

PRINT Moeaedededeacdededcdcdedededeaede dea sede! 

PRINT "WOULD YOU LIKE TO GET A HARD COPY?(Y/N)" 
QS=INKEYS: IF Q$="' " THEN 9280 

IF QS="N" OR QS$="n" THEN GOTO 10080 

IF QS="Y" OR QS="y'' THEN GOTO 9360 


GOTO 9280 
REM wevedevedeuvdedededeuesededededevede doves de de ves donk deve deve deve deve dev deve 
REM * PLOT THE LAST WAVEFORM * 


REM wrevedenkdedcdedevededevedede deve deve deve deuce He exe ae ve eve de deve de deck He deve ve 
VEWS="': HARD: SOUR FACT ,WMEM4" 


GOSUB 5840 
PRINT "PLOT ? (Y/N)" 
Q$S=INKEYS: IF QS=""_ "' THEN 9480 


IF QS="N" OR QS="'n'"' THEN GOTO 9600 

IF QS="Y" OR QS="y'"' THEN GOTO 9560 

GOTO 9480 

VEWS="': PLOT?;"’ 

GOSUB 5840 

GOTO 9780 

REM wevesedodkdvedodkvedede deve devede deve de veve sede de de dea deve He deve He se Fee ve He ae He He 


REM * PRINTING * 
REM  trvvededevededevedededeve vedo ve teveve donk ve see desk deve Te deve dee ve deve ak seve He ve 
VEWS=": HARDCOPY: PRINTER DEFAULT" 

GOSUB 5840 

VEWS="': HARDCOPY: PAGE AUTOMATIC” 

GOSUB 5840 

VEWS="': PRINT?; " 

GOSUB 5840 

BOADS="'GPIBO"' 

CALL IBFIND( BOADS , BOA%) 

CALL IBSIC( BOA%) 

VEWS="'7Q! '! 

CALL IBCMD (BOA%, VEWS ) 

V%=1 

CALL IBGTS (BOA%,V%) 

PRINT [lvecededkoedevedcdededededkdedsdededkdeukdededeaedededededk dese de desea a dete 
PRINT "PLEASE WAIT UNTIL THE PLOTTER HAD FINISHED" 
PRINT "PLEASE RELEASE THE DPO TO LOCAL(Y)"” 


QS=INKEYS: IF QS=" '' THEN 9280 
IF QS="Y" OR Q$="y" THEN GOTO 10040 
GOTO 9980 


CALL IBLOC(DSO0%) 
CALL IBLOC(BOA%) 


RETURN 
REM Jedeveteekioveleloitedchkkkcktlohketckkikhkidchkkhkkkstak 


REM SUBROUTINE READ DATA TO STRING REF. GURU D-47 
CALL IBRD(DS0%, RDS) 

RETURN 

REM wvedededevevedededevede devote dete de tedesedek KIT RR TRI I I 

REM * PROGRAM TO ARRANGE FORMAT FILE * 

REM * SCALE THE DATA ve 

REM * * 

REM * ¥e 


74 


10280 
10300 
10320 
10340 
10360 
10380 
10400 
10420 
10440 
10460 
10480 
10500 
OS20 
10540 
10560 
10580 
10600 
10620 
10640 
10660 
10680 
10700 
10720 
10740 
10760 
10780 
10800 
10820 
10840 
10860 
10880 
10885 
10890 
10892 
10896 
10897 
10898 
10899 
10900 
10901 
10902 
10903 
10904 
10910 
HOIZO 
10940 
10960 
10980 
11000 
11020 
11040 
11060 
11080 
11100 
PTZ 0 
11140 


REM Jevevedeve weve sede dees se veces teak vee vies ake ake ve ve ve se ve Fe 0 ve Fe 7 He Fe Ne He He Me 


OPEN "C: HEADER" FOR INPUT AS#1 
INPUT#1 ,A%,B%,C%,D%,E,F,G%,H,1,J,K 


PRINT "FORMAT "; A% 
PRINT "TYPE ";B% 
PRINT "POINT "30% 
PRINT ‘COUNT "; D% 
PRINT ''XINCRE ";E 
PRINT "XOR ay 
PRINT "XREF "3 G% 
PRINT "YING ";H 
PRINT "YOR eal 84 
PRINT "YREF ";J 
PRINT "RANGE ":K 
CLOSE#1 

REM DSES$="C:c" 
FOR K=1 TO ND 


IF K<10 THEN DSFILES=DSEWS+RIGHTS(STRS(K),1)+".wfm'" 
IF K>9 THEN DSFILES=DSEWS+RIGHTS(STRS(K),2)+'.wfm' 
IF K>99 THEN DSFILES=DSEWS+RIGHTS(STRS(K),3)+". wfm'' 
OPEN DSFILE$ FOR INPUT AS #1 


REM 
REM 
REM 
REM 
REM 
REM 
FOR I=1 TO C% 
INPUT #1,WFM 
WCI)= (CWFM-J)*H+IY) 
NEXT I 
SUMM=0. 0 
FOR I=1 TO C% 
SUMM=SUMM+W(T) 
NEXT I 
REM tetevedededscededevedede de cece rede Lee TOK OTT TE TERETE RTS TRICE 
iS dy as REMOVE MEAN ¥ 
REM wevevevedevedededsdedevede te se dete Fete Fe tent Fede Fe THe TCT REEL E 


SMEAN=SUMM/C% 

PRINT "MEAN ="; SMEAN 
HORsi—1 TOVCz 

RRR( I )=W( 1) -SMEAN 

NEXT I 

CLOSE#1 

DSFILES=DSE$+''. WFM" 

OPEN "A" ,#1,DSFILES 

PRINT#1,T$ 

PRINT#1,WS 

PRINT#1,D$ 

PRINT#1 , DUMMYS 

PRINT}##1 , D% 

PRINT#1,USING "sdf. #5 F1 

PRINT#1,ND 

PRINT#1,P1 

FOR I=1 TO C% 


PRINT#1, USING #4. {HEHEHE "; RRR(I) 


fe. 


11160 
11150 
11200 
PZ ZO 
Mbeya 
11240 
11260 
11280 
11300 
ali 3) 
11340 
11360 
11380 
11400 
11420 
11440 
11460 
11480 
11500 
11520 
11540 
IL 1bsyey 6. 
11580 
11600 
11620 
11640 
11660 
11680 
11700 
117 20 
11740 
bl 60 
WEF 3310. 
11800 
11820 
11840 
11860 
11880 
11900 
bo 28 
11940 
11969 
abe ren®, 
12000 
2020 
12040 
12060 
12080 
12100 
12120 
12140 
VA No) 
IZ 180 
12290 
12220 
12240 


NEXT I 
REM 

CLOSE# 
NEXT K 


1 


BEEP: BEEP 


RETURN 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 
REM 


Jv 


TSS=SPACES(64) 
DUMMYS=SPACES$( 64) 

KEY OFF: SCREEN 0,0,0 

OPEN "DATA. DAT' FOR INPUT AS #1 
INPUT#1,P1: INPUT#1,N: INPUT#1,ND: INPUT#1,T$: INPUT#1,D$: INPUT#1,T2 
INPUT#1 ,WS: INPUT#1,F1: INPUT#1,S1: INPUT#1,DSE$: INPUT##1, DUMMYS 
INPUT#1, DELAY1: CLOSE#1 
COLOR 7 ie oc us 


X=4: Y= 


16 


PRINT''PROGRAM DATA ACQUISTION " 


Mesdededesede sede dese dese se ve sede ve Teves Fete Teste ve se Ve He Fe Fe Fev He Fe 7 He Hee Ve Fete He tere Tee 
Ea = NUMBER OF POINTS 
N = NUMBER OF SUBAVERAGES 
ND = NUMBER OF DATA BLOCKS 
is = TARGET IDENTIFICATION 
DS = DATE OF MEASUREMENT 
he = DPO SAMPLING INTERVAL 
WS = WAVEFORM TYPE 
Fl = DPO TIME WINDOW 
Sl = DPO MAXIMUM SCALE 
DSES = OUTPUT DATA FILE NAME 
DUMMY$ = DUMMY § 


Jorkvoeveskosdeovotrscveseseseseskdedesedeskeakesesede sk deskscskseaksesdevedeakcseakseskesescae de ae 


FILE NAME FOR DATA OUTPUT...!... 


e4eoeeeeeeereeeoeoeoee#ee#eeeeeeee¢ 


NUMBER OF SUB AVERAGES.......... 


TARGET ID. (S$)... 30:0 0=oee t 


«ee @¢e8¢@8# e @ 0@ @© @ @ @ @ @ @ © @ @ @© @ @© @ @ 


WAVEFROM TYPE ! 


DPO SAMPLING INTERVAL.......... 
RUN THE PROGRAM..........e0eeee * 
DPO MAXIMUM VERT MV............ 
AUTO TIME@SET. ame. coe ee 


LOCATE X,Y+15: PRINT''PRESS THE FOLLOWING KEY" 
LOCATE X+2,Y: PRINT " < F > 

LOCATE X+3,Y: PRINT ' < P > NUMBER OF POINTS 
LOCATE X+4.YoRRING  < A> 

LOCATE X+5,Y¥: PRINT " < B > NUMBER OF DATA BLOCKS 
LOCATE X+6,Y: PRINT " < T > 

LOCATE X+7,Y: PRINT " < D> DATE ($) 

LOCATE X+8,Y: PRINT "' < W > 

LOCATE X+9,Y: PRINT " < K > DPO TIME WINDOW 
LOCATE X+10,Y: PRINT " < § > 

LOCATE X+11,Y: PRINT " < R > 

LOCATE X+12,Y:PRINT " < U > 

LOCATE X+13,Y: PRINT " < X > 

LOCATE X+14,Y: PRINT ' < G > DELAY TIME 
QS=INKEYS: IF QS=" " THEN 12000 

IF QS="G" OR QS$="g" THEN 12900 

IF Qs="F" OR QS="£"" THEN 12300 

IF QS="P" OR QS="p" THEN 12340 

IF QS="A" OR QS="a'"' THEN 12560 

IF QS="B" OR QS="'b" THEN 12600 

IF QS="'T" OR QS="t" THEN 12640 

IF QS="D" OR QS="d" THEN 12680 

IF QS="W" OR QS="'w" THEN 12700 

IF QS="S" OR QS="'s" THEN 12740 

IF QS="K" OR QS="k" THEN 12780 

IF QS="R" OR QS="r" THEN 4040 

IF Qs="'U" OR QS="u" THEN 12820 
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12260 
12280 
12300 
12320 
12340 
12360 
12380 
12400 
12420 
12440 
12460 
12480 
12500 
bz 20 
12540 
12560 
12580 
12600 
12620 
12640 
12660 
12680 
12700 
27 20 
12740 
oO) 
2760 
12800 
£2820 
12840 
12860 
12880 
12900 
E2920 
12940 
12960 
12980 
13000 
13020 
13040 
13060 


IF QS="X" OR QS="x" THEN 12860 

GOTO 11720 

PRINT" INPUT FILE NAME($) EX..C:F" 
INPUT DSES$: GOTO 12940 

CLS: LOCATE 4,8: PRINT''NUMBER OF POINT "' 
PRINT ''A)... POINT. "' 
PRINT "B).... POINT. "' 
PRINT ©). POINT." 
PRINT "D).. POINT." 
QS=INKEYS: 
16 @6= Ar 
5 Qs="B"' 
11 eS ou 
18 j= 0. 
GOTO 12940 
PRINT INPUT NUMBER OF SUB-AVERAGE" 

INPUT N: GOTO 12940 

PRINT 'INPUT NUMBER OF DATA BLOCK" 

INPUT ND :GOTO 12940 

PRINT INPUT TARGET ID(S$)" 

INPUT TS:GOTO 12940 

PRINT SET DATE": DS=DATES: GOTO 12940 

PRINT INPUT WAVEFROM TYPE(S)" 

INPUT WS: GOTO 12940 

PRINT" INPUT DPO SAMPLING INTERVAL" 

INPUT T2:GOTO 12940 

PRINT"DPO TIME WINDOW a 

INPUT F1:GOTO 12940 

PRINT" INPUT DPO VERTICAL SCALE" 

INPUT S$1:GOTO 12940 

PRINT'AUTO TIME" 

DUMMYS=TIMES: GOTO 12940 

PRINT INPUT TIME DELAY" 

INPUT DELAY1: GOTO 12940 

OPEN "DATA. DAT" FOR OUTPUT AS#1 

PRINT#1,P1: PRINT#1,N: PRINT#1,ND: PRINT#1,TS: PRINT#1,D$ 
PRINT#1,T2: PRINT#1,WS: PRINT#1,USING "sy. dH"; F1: PRINT#1,S1 
PRINT#1,DSES: PRINT##1 , DUMMYS: PRINT#1,DELAY1 

CLOSE#1 : CLS 

REM 

GOTO 11700 


OR QS="'b" THEN P1=512 
OR Q$="'c" THEN P1=256 
OR QS="'d"' THEN P1=128 
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APPENDIX C. CHARACTERISTICS OF LOW LOSS CABLE 


Table 21. CHARACTERISTICS OF PORT RETURN CABLES 


















Dielectric withstand LOkV 


Table 22, INSERTION LOSS 
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